AD-A031  057 


OHIO  STATE  UNI V RESEARCH  FOUNDATION  COLUMBUS  F/G  20/8 

NUCLEAR  MAGNETIC  RESONANCE  INVESTIGATION  OF  ELECTRON  DAMAGE  IN  — ETC(U) 


1 


APR  76  L C BROWN  F33615-74-C-2058 

UNCLASSIFIED  QSURF-3909 AFAPL-TR-76-36  , NL 


1 0F2 
AO 

AO  3 105  7 

a 

& 

- 

I—— 

r— 

1 

\r 

' 

MICROCOPY  RESOLUTION  TEST  CHARI 

NATIONAL  BURtAU  OF  STANDARDS -1963-4 


AFAPLTR-76-36 


NUCLEAR  MAGNETIC  RESONANCE 
INVESTIGATION  OF  ELECTRON  DAMAGE  IN 
□ DOPED  SILICON 


THE  OHIO  STATE  UNIVERSITY 
RESEARCH  FOUNDATION 
isn  KINNEAR  ROAD 
COLUMBUS,  OHIO  J>3212 


APRIL  1976 


TECHNICAL  REPORT  AFAPL-TR-76-36 

FINAL  REPORT  FOR  PERIOD  16  MAY  1974  - 15  MARCH  1976 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  AERO-PROPULSION  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT  PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a definitely  related  government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in  any  way  supplied  the 
said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

This  final  report  was  submitted  by  The  Ohio  State  University  Research 
Foundation,  under  Contract  F33615-74-C-2058.  The  effort  was  sponsored 
by  the  Air  Force  Aero  Propulsion  Laboratory,  Air  Force  Systems  Command, 
Wright-Patterson  AFB,  Ohio,  with  W.  Patrick  Rahilly/AFAPL/POE-2/255-6237 
as  Project  Engineer  In  Charge.  Prof.  L.  Carlton  Brown  of  The  Ohio  State 
University  was  technically  responsible  for  the  work. 

This  report  has  been  reviewed  by  the  Information  Office,  (ASD/OIP) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS). 

At  NTIS,  it  will  be  available  to  the  general  public,  including  foreign 
nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publi- 
cation. 


Name  and  Grade 

Project  Engineer/Scientist 

FOR  THE  COMMANDER 


Name  and  Title 


N 

S 

m 

i 


Copies  of  this  report  should  not  be  returned  unless  return  is  re- 
quired by  security  considerations,  contractual  obligations,  or  notice  on 
a specific  document. 


AlW  FOHCF  - 7?  Au4f.T  7 f>  - 150 


Unclassified 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (When  Data  Entered) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
3.  RECIPIENT’S  CATALOG  NUMBER 


/^.REPORT  DOCUMENTATION  PAGE 


4.  TITLE  (and  Subtitle) 


„ NUCLEAR  J1AGNETIC  RESONANCE  JNVESTIGATION  OF 
"'ELECTRON  DAMAGE  IN  Li-DOPED 'SILICON  . 


Final 


NUMBER(S) 


7 AUTHORS 


L . Carlton, 


irown 


10  program  element,  project,  t ask 


9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


The  Ohio  State  University  , 

Research  Foundation 

1314  Kinnear  Road,  Columbus,  Ohio  4-321 2 


CONTROLLING  OFFICE  NAME  ANO  ADDRESS  ^ 

Department  of  the  Air  Force  ' j j 

Air  Force  Aero  Propulsion  Laboratory  (AFSC)  ' — 
Wright-Patterson  Air  Force  Base,  Ohio  4-5433 

. MONITORING  AGENCY  NAMF.  ft  AODRESSf//  different  from  Controlling  Office) 


13.  NUMBER  OF  PAGES 

125 

15.  SECURITY  CLASS,  (of  this  report) 


15a  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited, 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  mnrher) 


nuclear  magnetic  resonance  (NMR) 
n-type  crucible-grown  Li-doped  silicon 
electron  irradiation  damage 
electron  paramagnetic  resonance  (EPR) 


lithium-doped  silicon 


20  ABSTRACT  (Continue  on  reverse  aide  if  necessary  and  identify  bv  bloc*  numoer ) 


FORM 
1 JAN  73 


Unclassified 

SECURITY  CLASSlF | CATION  OF  THIS  PAGE  fWhen  Data  Entered' 


EDITION  OF  1 NOV  65  IS  OBSOLETE 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE('W7»«n  Data  Entered) 


Block  #20  (Abstract) — Continued 

deep  levels.  In  upper  temperature  ranges,  where  the  conduction  electron 
concentration  is  the  dominant  factor,  the  relaxation  rate  follows  changes 
in  conduction  electron  concentration  resulting  from  changes  in  donor  impur- 
ity concentration.  At  somewhat  lower  temperatures,  where  the  conduction 
electron  spends  part  of  the  time  in  a trap  level,  the  relaxation  rate  shows 
a dependence  upon  the  NMR  operating  frequency.  The  relaxation  rate  is 
higher  at  lower  operating  frequencies,  indicating  that  the  electron  resi- 
dence time  (or  correlation  time)  is  longer  than  the  Larmor  period  for  one 
radian  of  nuclear  precession  for  the  magnetic  fields  employed.  At  even 
lower  temperatures,  below  4o  K or  50  K,  the  electrons  are  essentially  all 
localized  and  the  NMR  relaxation  rate  becomes  dominated  by  the  correlation 
time  of  the  trapped  paramagnetic  electron  in  one  of  its  spin  states,  which 
is  closely  related  to  the  relaxation  time,  Te,  of  the  electron  spin.  The 
NMR  relaxation  rate  increases  sharply  with  decreasing  temperature  between 
50  K and  20  K and  becomes  much  less  dependent  upon  NMR  operating  frequency, 
indicating  that  Te  is  smaller  than  the  Larmor  period.  At  about  30  K-40  K 
the  steady  state  NMR  magnetization,  Ms,  which  was  also  measured  in  the 
study,  was  found  to  deviate  from  Curie-law  behavior,  to  decrease  at  lower 
temperatures,  and  to  disappear  at  about  20  K in  all  except  the  lithium-rich 
10-17  sample.  In  the  latter  sample  Mg  decreased,  leveled  off,  and  then  be- 
gan to  increase  with  decreasing  temperature  at  about  10  K (before  irradia- 
tion). This  behavior,  along  with  results  obtained  from  EPR  and  Hall  mea- 
surements, also  a part  of  the  study,  supports  the  conclusion  that  lithium 
forms  pairs  or  clusters,  with  itself  or  with  other  impurities,  which  tend 
to  become  diamagnetic  at  low  temperatures  through  electron  exchange  inter- 
actions. / 

Changes  in  the  variation  of  the  NMR  relaxation  rate  vs  temperature 
after  electron  irradiation  (~  1 MeV  at  1016  cm-2  fluence)  and  again  after 
annealing,  as  well  as  changes  in  other  measured  parameters,  NMR  magnetiza- 
tion Ms,  EPR  magnetization  Me,  and  Hall  effect  related,  indicate  an  increase 
in  the  variety  of  trap  levels  after  irradiation  with  only  a slight  reduction 
in  variety  after  annealing.  In  most  samples  irradiation  appeared  to  produce 
more  deep  levels  than  shallow  (as  compared  with  phosphorus),  however  in  the 
lithium-rich  samples  more  shallow  levels  appeared  to  be  produced  than  deep 
levels. 
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SUMMARY 


SOLAR  CELL  MATERIALS 

Bulk  sample  measurements  have  been  made  on  n-type  crucible-grown 
silicon  containing  three  different  phosphorus  concentrations  with  10, 
1,  and  0.1  fi-cm  starting  resistivities.  Some  of  the  sample  wafers 
(20  x 10  x 0.25  mm)  also  contained  lithium  in  several  concentrations 


(of  -lO59,  1013,  101 


1019  cm-3), 


The  samples  studied  by  both 


NMR  and  EPR  techniques  were  10-15,  10-17,  1-0,  1-15 , 0.1-15  where  the 
first  number  refers  to  the  starting  resistivity  and  the  second  is  the 
exponent  of  the  Li  concentration.  EPR  investigations  have  been  extended 
to  other  samples. 


EXPERIMENTAL  PARAMETERS  CONTROLLED  AND  MEASURED 

In  the  NMR  investigation  the  29Si  Spin- lattice  relaxation  time,  Ti , 
and  the  steady  state  nuclear  radiofrequency  magnetization  magnitude,  Ms, 
were  measured  and  the  spin-spin  relaxation  time  T 2 (related  to  line 
width)  was  monitored  as  a function  of  NMR  operating  frequency  (field) 
and  sample  temperature,  T,  over  a range  from  8 to  310  K.  The  NMR  operat- 
ing frequency,  f,  was  set  at  approximately  6,  8.5,  and  10. 5 MHz  in  order 
to  obtain  the  effects  of  a variation  in  operating  frequency  (field). 

The  measurements  indicated  above  were  made  on  each  of  the  samples  in- 
vestigated, first  in  their  original  state,  then  again  after  irradiation 
with  ~1  MeV  electrons  to  a fluence  of  10 16  electrons/cm2.  Immediately 
after  irradiation  the  samples  were  stored  at  77  K to  prevent  Li  diffu- 
sion until  the  second  round  of  measurements  was  completed.  Finally  the 
measurements  were  repeated  after  annealing  the  samples  at  room  tempera- 
ture for  periods  of  12  to  60  hours . Each  NMR  sample  consisted  of  15  of 
the  20  x 10  x 0.25  mm  wafers  or  0.75  cm3.  NMR  measurements  of  Tx  were 
made  for  temperatures  to  400  K for  some  samples. 

In  the  EPR  investigation  a smaller  sample  4 x 10  x 0.25  mm  or 
0.01  cm3  was  used  and  measurements  were  made  before  and  after  1 MeV 
electron  irradiation  (fluence  of  1016  electrons /cm3)  and  after  corre- 
sponding annealing  at  room  temperature.  The  measurements  were  made  at 
X-band  microwave  frequencies  and  at  temperatures  ranging  from  6 to  310  K. 
The  quantities  measured  were  the  effective  g factor  (spectral  position), 
electronic  steady  state  magnetization  magnitude,  Mg,  line  width  (re- 
lated to  electronic  spin-spin  relaxation,  T2e)>  and  splitting  factor, 

AG,  in  the  cases  where  multiplet  structure  was  observed. 


Supporting  Hall  coefficient  and  resistivity  data  were  measured 
under  the  same  conditions  of  temperature,  irradiation,  and  annealing. 

A 20  x 10  x 0.25  mm  sample  size  was  used  with  five  contacts  of  gallium- 
indium  alloy  applied  with  ultrasonic  vibrations. 


EXPERIMENTAL  RESULTS  AND  FHYSICAL  MODELS  - NMR 


In  measuring  the  r3Si  spin-lattice  relaxation  time,  Tj. , a saturating 
operation  is  first  performed  on  the  sample  in  order  to  reduce  the  magne- 
tization M,  due  to  the  ?9Si  nuclear  moments,  to  zero;  i.e.,  the  spin 
populations  are  equal  for  the  m = \ and  m = spin  states . Then  the 
NMR  field  and  sample  temperature  are  held  constant  for  a period  of  time 
t during  which  the  magnetization  M is  found  to  grow  exponentially  accord- 
ing to  the  relation  M(t)  = Ms[l  - exp(-t/Ti)].  The  state  of  M(t)  is 
tested  at  time  t by  means  of  a quasi-adiabatic  fast-passage  through  the 
NMR  resonance,  sensing  the  dispersion  mode.  The  process  is  repeated  for 
several  values  of  t.  The  quantity  l/Ti  then  represents  the  rate  of 
recovery  or  growth  and  Ms  represents  the  steady-state  or  maximum  value 
of  M. 


In  this  investigation  the  rate  of  recovery  l/Tx  is  controlled 
entirely  by  either  (l)  conduction  electrons  or  (2)  paramagnetic  elec- 
trons such  as  single  electrons  occupying  a shallow  or  deep  trap.  There 
are  two  main  factors  involved  in  the  influence  of  the  electrons  on  l/TL . 
One  of  these,  the  strength  of  the  electron-nuclear  interaction,  can  be 
thought  of  as  the  strength  of  the  local  magnetic  field,  due  to  the 
electron,  at  the  site  of  a >9Si  nucleus.  This  field  is  called  hj^t). 

Of  course,  different  nuclei  will  "see"  different  values  so  that  hp 
represents  a mean  of  a statistical  variable  with  some  distribution  about 
the  mean.  The  second  factor  in  the  influence  of  the  electron  on  l/Tj. 
is  the  time  variation  of  hp(t)  or  specifically  the  component  in  the  fre- 
quency spectrum  (power  spectrum)  of  hp(t)  which  falls  at  the  NMR  operat- 
ing frequency,  cj  = 2nf. 


The  power  spectrum  of  the  electronic  local  field  is  related  to  the 
auto-correlation  function  of  hp(t)  through  the  quantity  t called  the 
auto -correlation  time  (or  correlation  time).  A number  of  theories  have 
been  treated  relating  to  this  problem,  of  which  the  most  reasonable 
ones  find  that  the  power  spectrum  and  therefore  l/Tv  are  proportional 
to  hpt(l  + w t ) . The  quantity  t can  then  be  interpreted  as  a measure 
of  the  lifetime  of  an  electron  in  a localized  state  or  in  a spin  state. 
The  frequency  w is  then  taken  to  be  the  NMR  operating  frequency. 


The  above  relation  appears  in  two  different  temperature  ranges  in 
the  analysis  of  the  experimental  results.  In  the  range  from  310  K down 
to  about  60  K,  t is  interpreted  as  the  lifetime  of  an  electron  in  a 
trap.  At  temperatures  below  about  40  K,  t is  related  to  the  relaxation 
time  of  the  electron;  i.e.,  the  lifetime  in  a spin  state. 


Upper  Temperature  Region 


In  the  upper  temperature  range  the  lifetime  of  an  electron  in  a 
spin  state,  Te,  is  very  short,  as  is  seen  by  the  very  wide  EPR  linewidth 
observed  at  310  K (Te  less  than  lO-^  s).  Under  these  conditions  the 
^i  nucleus  "sees"  the  average  value  of  the  electron  magnetic  moment, 
averaged  over  m \ and  m -'  spin  states,  however,  the  average  value 


is  not  zero,  but  rather  u = ne(beHoA-T)  where  the  fraction  in  parentheses 
comes  from  the  Boltzmann  factor . The  correlation  time  t in  this  range 
then  represents  the  longer  time  that  an  electron  sits  in  a trap  and 
appears  (to  the  ^Si)  to  have  a moment  p (less  than  pe)  and  a local  field 
h^  proportional  to  p.  In  the  plots  of  l/Tx  vs.  T (see  Pigs.  5-19). 
samples  0.1-15  and  10-17,  dominated  by  phorphorus  and  lithium,  respec- 
tively, the  measurements  are  seen  to  be  frequency  independent  (converge) 
at  room  temperature,  to  drop  to  lower  rates  as  the  temperature  is 
lowered,  and  to  become  increasingly  frequency  dependent.  This  is  an 
indication  that  wt  beecmes  increasingly  greater  than  unity  or  t > l/w. 

For  the  NMR  frequencies  of  6,  8.5,  and  10.5  MHz  l/w  has  the  values  of 
28,  20,  and  16  ns,  respectively,  thus  providing  lower  limits  to  the 
localized  electron  lifetimes.  Below  about  60  K a different  mechanism 
becomes  dominant  as  discussed  below. 

A comparison  of  the  plots  for  0.1-15  and  10-17  shows  that  the  10-17 
plot  has  a smaller  slope  (in  the  60  to  310  K region),  is  shifted  slightly 
toward  lower  temperatures,  and  shows  a slight  frequency  dependence  at 
room  temperature . This  is  interpreted  to  be  the  result  of  a smaller 
average  trap  binding  energy  for  Li  (lower  temperature)  and  a broader 
distribution  of  trap  depths  about  the  average  (smaller  slope)  including 
some  larger  trap  energies  (room  temperature  frequency  dependence),  thus 
supporting  the  contention  that  Li  in  oxygen-rich  silicon  forms  several 
different  types  of  donor  sites. 

In  the  samples  of  lower  donor  concentrations  (10-15,  1-15,  1-0) 
the  room  temperature  frequency  dependence  is  proportionally  larger  as 
compared  with  the  overall  smaller  rate  l/Tj_  (longer  TL)  indicating  that 
deeper  traps  than  phosphorus  are  competing  on  a roughly  equal  basis  in 
contributing  to  ?9Si  relaxation.  A projection  of  the  curves  to  higher 
temperatures  suggests  that  perhaps  600  K would  be  required  to  reach 
frequency  independence,  implying  that  traps  perhaps  twice  as  deep  as 
phosphorus  are  present  in  roughly  equal  concentrations.  Since  the 
slopes  are  smaller,  a distribution  of  trap  energies  is  indicated. 

The  plot  of  l/Tj.  vs.  T for  other  samples,  before  and  after  irra- 
diation and  annealing,  have  been  studied  in  a similar  manner  and  the 
indications  in  the  upper  temperature  range  (60  to  310  K)  are  that,  gen- 
erally, irradiation  causes  a small  decrease  in  temperature-independent 
room  temperature  relaxation  rate,  corresponding  to  a decrease  in  con- 
duction electron  density,  followed  by  a smaller  increase  after  anneal- 
ing. An  increase  in  frequency  dependence  and  a smaller  slope  is  evident 
after  irradiation  and  this  effect  does  not  recover  after  annealing. 

This  indicates  that,  generally,  there  is  a broader  distribution  of  trap 
energies  after  irradiation  which  persists  after  annealing. 

Low  Temperature  Region 

In  the  lower  temperature  region  (20  to  40  K)  there  is  a sharp  in- 
crease in  rate  l/T^ , with  decreasing  temperature,  observed  in  all  samples. 
Since  the  electrons  have,  at  this  temperature,  long  lifetimes  as 


v 


localized  electrons,  the  power  spectrum  in  this  range  results  from  the 
electron  changing  spin  states  or  from  ground  state  to  excited  state  in 
the  trap  potential.  A number  of  investigators  have  studied  the  electron 
relaxation  rate  l/Te  in  this  region  and  have  found  that  different 
processes  have  temperature  dependence  proportional  to  T or  T9  or 
exp  A/kT.  These  studies  all  indicate  that  there  is  a very  rapid  in- 
crease in  Te  with  decreasing  temperature,  and,  if  this  is  taken  to  be 
the  correlation  time  t in  the  expression  hp-r/^l  + u)pt2)  for  the  power 
spectrum  of  the  local  field,  a useful  model  results  for  the  interpreta- 
tion of  l/Tx  vs.  T in  this  region.  The  position  of  the  characteristic 
slope  in  this  region  is  shifted  toward  lower  temperatures  in  the  case 
of  smaller  trap  energies  (as  in  10-17)  and  shows  a smaller  slope  in  the 
case  of  a distribution  of  trap  energies.  A study  of  the  plots  of  l/TL 
vs.  T in  this  range  again  shows  an  increase  in  average  trap  energy  and 
a broader  distribution  of  energies  after  irradiation  which  does  not 
recover  completely  after  annealing. 

NMR  Magnetization  vs . Temperature 

Elementary  theory  indicates  that  the  NMR  magnetization  (steady 
state),  Ms,  should  increase  with  decreasing  temperature  in  proportion 
to  l/T  (Curie's  Law).  The  plot  of  Ms  vs.  T (Figs.  4,  20)  shows  Curie's 
Law  behavior  down  to  temperatures  of  approximately  40  K below  which  a 
sharp  drop  in  Ms  occurs,  with  different  samples  dropping  off  at  different 
temperatures . A comparison  of  the  10-15,  l-15j  and  0.1-15  reveals  that 
higher  P concentrations  drop  at  higher  temperatures.  In  contrast,  a 
comparison  of  10-17  vs.  (10-15,  1-15>  0.1-15)  vs.  1-0  reveals  that  higher 
Li  concentrations  drop  off  at  lower  temperatures.  This  indicates  (as 
expected)  that  larger  phosphorus  concentrations  correspond  to  higher 
concentrations  of  (unpaired)  paramagnetic  electrons  at  low  temperatures 
which  serve  to  remove  contributions  to  Mg  from  the  observed  NMR  line 
when  seme  of  the  electrons  have  t (or  Tej  sufficiently  long  (in  the 
millisecond  range)  to  produce  locally  inhomogeneous  fields. 

The  surprising  result  in  the  contrasting  behavior  of  Ms  vs.  Li 
concentration  is  taken  as  evidence  that  Li,  being  mobile,  tends  to  form 
Lis  pairs  (and  also  Li  P pairs)  or  clusters  which,  at  low  temperatures, 
become  partially  diamagnetic  thus  reducing  the  low  temperature  paramag- 
netic concentration  with  increasing  Li  concentration.  These  conclusion 
are  supported  by  evidence  obtained  from  EPR  results  and  from  Hall  effect 
and  mobility  studies.  This  effect  is  particularly  noticeable  in  the 
10-17  samples  where  Ms  (with  decreasing  temperature)  starts  to  decrease 
and  then  levels  off  and  starts  to  increase  at  10  K,  indicating  that 
many  of  the  paramagnetic  sites  at  20  K have  become  diamagnetic  at  10  K . 


The  behavior  of  Ms  vs  T with  varying  lithium  concentration,  as  con- 
trasted with  varying  phosphorus  concentration,  is  thus  attributed  to  the 
mobility  of  the  Li  ion  and  to  a tendency  of  the  Li  ion  to  form  a weak 
bond  with  other  donor  impurity  sites  (including  Li -Li,  Li-LiOV,  Li-P, 
Li-deep  trap,  and  others)  partially  due  to  electron  exchange  interactions 
as  discussed  in  the  EPR  section  in  some  detail  (see  Eq.  86  and  following). 
The  most  probable  situation  for  Li,  being  mobile,  is  to  minimize  the 


I 


total  energy  and  thus  to  maximize  the  value  of  the  exchange  energy 
parameter  J by  adjusting  the  distance  from  neighboring  donors  (pairing 
or  clustering).  This  situation  also  serves  to  maximize  the  population 
of  the  lowest  level  singlet  state  (S  = 0).  Since  this  state  is  diamag- 
netic, the  paramagnetic  electron  population  at  lower  temperatures  is  re- 
duced below  that  expected  in  the  case  of  a pure  system  of  immobile  ions, 
such  as  phosphorus,  which  cannot  make  such  adjustments  in  position.  All 
of  the  donor  sites  do  not  become  diamagnetic,  of  course,  and  those  that 
remain  paramagnetic , whether  isolated  sites  or  clusters  with  nonzero 
spin,  eventually  attain  correlation  times  and  effective  magnetic  moments 
of  sufficient  magnitude  to  displace  neighboring  29Si  spins  from  the  ob- 
served NMR  signal  and  thus  cause  a reduction  of  Ms,  however,  at  lower 
temperatures  than  in  lithium- lean  samples. 


EXPERIMENTAL  RESULTS  AND  PHYSICAL  MODELS  - EPR 

At  room  temperature  the  EPR  line  is  very  weak  and  broad,  if  seen 
at  all,  indicating  a very  short  lifetime  of  the  electron  in  a paramag- 
netic state  (one  ns  or  less).  At  about  100  K the  line  has  narrowed  to 
about  3 G (Tpe  about  10  ns)  and  the  g-factor  can  be  measured  to  better 
than  one  part  in  103.  In  P-rich  samples,  g corresponds  to  reported 
values  for  phosphorus  site  electrons.  In  Li-rich  samples,  two  different 
g factors  were  observed  and  were  in  agreement  with  g for  LiO  and  LiOV 
type  sites. 

At  temperatures  lower  than  30  K the  P-site  line  was  observed  to 
split  into  a doublet  (aG  = 42  G)  and  in  P-rich  samples  weaker  triplet 
and  quadruplet  structure  was  observed.  This  splitting  behavior  results 
from  interaction  of  the  electron  with  the  phosphorus  nucleus  (s  = -y)  in 
the  case  of  the  42  G doublet.  The  triplet  and  quadruplet  result  from 
an  exchange  interaction  between  two  or  three  adjacent  phosphorus  sites 
and  appear  when  kT  has  become  smaller  than  the  exchange  energy  J.  The 
two  electrons  can  assume  a triplet  state  (S  = l)  and  remain  paramagnetic, 
or  a singlet  state  (S  = 0)  and  become  diamagnetic.  In  the  latter  case 
the  contribution  is  lost  from  EPR  magnetization  and  also  NMR  magnetiza- 
tion is  increased  and  NMR  relaxation  rate  is  decreased.  In  the  triplet 
state  the  two  phosphorus  nuclear  moments  combine  to  produce  the  three- 
component  line  in  the  EFR  spectrum.  The  lithium  lines  do  not  show  a 
nuclear  hyperfine  splitting  even  though  the  dominant  7Li  nucleus  has  a 
spin  of  3/2.  This  results  from  the  lower  symmetry  of  the  interstitial 
Li  position  and  the  strong  p-state  character  of  the  donor  electron  as 
contrasted  with  the  strong  s-state  character  of  the  phosphorus  donor 
electron.  In  the  p-state  the  electron  density  at  the  nucleus  is  small 
and  thus  also  the  hyperfine  interaction. 

The  pairing  of  Li  (or  clustering)  is  revealed,  however,  in  the 
behavior  of  the  10-17  sample  upon  irradiation . An  increase  in  electron 
magnetization  of  an  order  of  magnitude  was  observed  at  25  K after  irra- 
diation, which  then  reduced  nearly  to  the  original  level  after  24  hr 
room  temperature  annealing.  This  can  be  accounted  for  by  a break-up  of 
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lithium  clusters  by  the  irradiation,  followed  by  a relatively  rapid 
return  to  cluster  formation  upon  annealing.  This  model  is  also  sup- 
ported by  the  Hall  measurements  and  by  the  NMR  Ms  vs.  T data. 

In  contrast  with  Curie  Law  behavior,  all  of  the  samples  generally 
show  a decrease  in  electronic  steady-state  rf  magnetization.  Mg,  with 
decreasing  temperatures  in  the  region  below  about  25  K.  Several  mecha- 
nisms can  contribute  to  this  effect:  (l)  change  in  neighboring  pairs 

from  paramagnetic  to  diamagnetic  states  when  kT  becomes  less  than  the 
pair  exchange  energy  J,  (2)  increase  in  local  field  inhomogeneity  when 
the  electronic  relaxation  time  becomes  sufficiently  long,  or  (3)  partial 
saturation  of  the  spin  system.  Tests  were  made  to  see  if  partial  satu- 
ration could  be  present  but  the  indications  were  that  this  was  not 
responsible  for  the  observed  non-Curie  behavior.  Evidence  from  EPR, 

NMR,  and  Hall  measurements  all  tend  to  support  the  first  two  of  the 
above  mechanisms.  Thus,  if  two  (or  more)  donor  sites  are  neighbors 
(30-40  A or  less),  whether  they  are  P-P,  Li-Li,  P-Li,  Li-LiOV,  or  other, 
at  low  temperatures  the  electrons  localized  on  these  sites  form  exchange- 
energy  pairs,  some  of  which  are  in  the  diamagnetic  lowest  energy  state 
and  do  not  contribute  to  the  EPR  magnetization  Me.  In  addition,  the 
remaining  paramagnetic  electrons  achieve  long  relaxation  times  suffi- 
cient to  increase  the  local  field  inhomogeneity  and  remove  contributions 
from  the  observed  EPR  linewidth. 

The  increase  in  low- temperature  EPR  magnetization,  Me,  after 
electron  irradiation  is  interpreted  to  indicate  that  these  pairs  (or 
clusters)  are  scattered  or  broken-up  by  the  bombardment,  and  further 
that  the  clusters  tend  to  reform  after  annealing,  again  reducing  Mg. 


HALL  EFFECT  AND  RESISTIVITY  MEASUREMENTS 

Conduction  electron  density,  n;  resistivity,  p;  and  mobility,  p, 
were  obtained  from  Hall  effect  and  resistance  measurements  as  a function 
of  temperature  and  then  redetermined  after  electron  irradiation  and 
after  annealing  (see  Figures  22-25). 

Generally,  the  conduction  electron  density  showed  a decrease  after 
irradiation,  indicating  a loss  of  conduction  electrons  to  deep  traps, 
except  in  the  case  of  the  10-17  sample  (lithium-rich)  where  the  conduc- 
tion electron  density  showed  an  increase  after  irradiation  as  well  as 
after  annealing.  This  is  viewed  as  supporting  evidence  for  lithium 
pairing  or  clustering.  After  irradiation,  the  deeper  trap  defects  are 
outnumbered  in  the  lithium-rich  sample  by  the  scattered  Li  interstitial 
shallow  donors  thus  increasing  the  conduction  electron  density.  Upon 
annealing,  the  further  increase  in  conduction  electron  density  indicates 
that  the  lithium  heals  the  deep  trap  defects  more  .rapidly  than  it  re- 
turns to  clusters. 

The  mobility,  generally,  shows  an  increase  with  decreasing  tempera- 
ture, however  the  greatest  rate  of  increase  occurs  after  irradiation 
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(before  annealing),  especially  in  the  lithium-bearing  samples.  Since 
the  low-temperature  mobility  reflects  the  contribution  from  shallow- 


trap  electrons  (the  only  ones  still  ionized),  the  greater  low-temperature 
mobility  after  irradiation  apparently  results  from  very  shallow  traps, 
again  supporting  the  increase  in  isolated  Li  interstitial  donors  upon 
irradiation . 
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INTRODUCTION 


GENERAL 

-'^This  contract  effort  represents  an  experimental  and  analytical 
investigation  of  electron  irradiation-damage  and  lithium-annealing 
processes  in  bulk  n-type  crucible-grown  silicon  samples.  Objectives  of 
the  program  were  to  obtain  experimental  data  from  nuclear  magnetic 
resonance  (NMR),  electron  paramagnetic  resonance  (EFR),  and  Hall  effect 
measurements,  and  to  correlate  these  data  with  the  assistance  of  appro- 
priate accepted  theories,  all  with  the  goal  of  increasing  the  knowledge 
of  the  radiation-damage  and  lithium-annealing  processes.,— — 
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TECHNICAL  APPROACH 

In  the  NMR  investigation,  the  29Si  nuclear  magnetic  moments 
(spin  - -g)  were  used  as  effective  probes  to  reveal  energy  transfer  and 
other  processes  in  the  bulk  samples.  The  spin-lattice  relaxation  time, 
TL , and  the  equilibrium  nuclear  magnetization,  Ms,  were  measured  as  a 
function  of  temperature  from  8 to  310  K (to  400  K in  selected  cases). 
The  linewidth,  l/T2,  and  spectral  position  (Knight  shift)  were  also 
monitored.  These  measurements  were  made  at  three  different  NMR  operat- 
ing frequencies,  w,  before  electron  irradiation,  after  irradiation,  and 
after  annealing.  The  samples  were  irradiated  with  ~1  MeV  electrons  at 
room  temperature  to  a fluence  of  1016  cm*2  and  immediately  stored  in 
liquid  nitrogen  to  inhibit  lithium  diffusion  until  the  second  round  of 
measurements  was  completed.  The  samples  were  then  annealed  at  room 
temperature  for  periods  varying  from  12  to  60  hours  and  reexamined. 
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Three  different  phosphorus  concentrations  were  used,  corresponding 
to  resistivities  of  0.1,  1,  and  10  ft- cm.  The  lithium  concentrations 
used  were  0,  1015,  and  101"7  cm-3.  The  samples  studied  were  labeled 
10-17,  10-15,  1-15,  0.1-15,  and  1-0  where  the  first  number  refers  to 
the  phosphorus -related  resistivity  and  the  second  to  the  exponent  of 
the  lithium  density.  Each  NMR  sample  consisted  of  15  slabs,  20  x 10  x 
0.25  mm,  insulated  from  one  another  by  teflon  film,  to  provide  a silicon 
volume  of  0.75  cm3. 

In  the  EPR  investigation,  the  same  sample  categories  were  studied 
under  the  same  conditions  of  irradiation  and  annealing  and  over  the  same 
range  of  temperatures.  The  measured  quantities  were  rf  electronic 
equilibrium  magnetization,  Me;  spectral  position,  g;  linewidth,  l/T2e; 
and  splitting  factor  AG,  in  the  case  of  multiplet  structure.  The  sample 
size  was  10  x 4 x 0.25  mm  or  0.01  cm  . 

The  Hall  measurements  included  the  determination  of  conduction 
electron  density,  resistivity,  and  mobility  for  the  same  sample  cate- 
gories, the  same  conditions  of  irradiation  and  annealing,  and  over  the 
same  range  of  temperatures . 
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NMR  INVESTIGATIONS 


GENERAL 

The  usefulness  of  NMR  in  studying  semiconductor  properties  and 
radiation  damage  in  silicon  results  from  a number  of  circumstances. 

The  active  magnetic  nucleus,  °9Si,  accounts  for  4.67 $ of  the 
naturally  occurring  silicon,  the  remaining,  92.28$  PeSi  and  3-05$  3°Si, 
have  zero  nuclear  spin  and  are  inert  from  the  standpoint  of  NMR.  Then, 
there  is  a relatively  dilute  and  stable  active  nuclear  moment  system 
numbering  about  2 x 10 11  in  the  0.75  cm'  sample  volume  used  in  the 
experiment.  This  provides  a reasonably  narrow  (less  than  0.4  G)  line- 
width  in  an  external  magnetic  field  of  ~ 10,000  G and  a useable  signal - 
to-noise  ratio  for  the  dispersion  mode. 

The  9Si  system  concentration  is,  of  course,  unaffected  by  changes 
in  field,  sample  temperature,  electron  irradiation,  annealing,  and 
lithium  concentration  and  is  essentially  unaffected  by  changes  in 
phosphorus  concentration. 

The  ’Si  nuclear  spin  of  l/2  ensures  that  direct  transfer  of 
nuclear  magnetic  energy  (Zeeman  interaction  energy  between  nuclear 
moments  and  external  magnetic  field)  to  the  crystal  lattice  is  virtually 
impossible,  however,  it  is  easy  to  transfer  energy  into  the  spin  system 
by  means  of  a radiofrequency  magnetic  field  whose  frequency  is  adjusted 
to  the  Larmor  precession  frequency  of  the  99Si  nuclei  in  a strong  con- 
stant external  magnetic  field.  The  time  scale  (spin-lattice  relaxation 
time,  TL ) for  energy  transfer  from  the  spin  system  to  the  lattice  in 
high-purity  silicon  at  low  temperatures  has  been  estimated  to  be  of  the 
order  of  weeks  or  more.  However,  if  electrons  of  a paramagnetic  charac- 
ter are  introduced,  such  as  non-degenerate  conduction  band  electrons, 
degenerate  conduction  band  electrons  at  finite  temperatures,  or  localized 
unpaired  electrons,  in  traps  shallow  or  deep,  the  spin-lattice  relaxation 
time  of  the  nuclear  spin  system  is  reduced  to  the  order  of  tens,  hun- 
dreds, or  thousands  of  seconds  and  these  times  can  be  measured  with 
reasonable  success  in  the  laboratory.  Thus,  the  rate  of  transfer, 
l/Tj. , of  energy  from  the  nuclear  spin  system  to  the  lattice,  in  the 
specific  samples  investigated,  for  all  practical  purposes,  is  entirely 
dependent  upon  the  state  of  the  donor  electron  system,  whether  in  the 
conduction  band  or  in  shallow  or  deep  traps . 

The  conduction  band  electrons  contribute  to  this  rate  in  a dif- 
ferent way  from  the  localized  electrons  as  functions  of  sample  tempera- 
ture and  Larmor  frequency  so  that,  by  controlling  these  laboratory 
parameters,  the  separate  contributions  can  be  sorted  out  and  analyzed 
in  terms  of  irradiation  and  annealing  effects. 

The  theoretical  structure  that  describes  nuclear  magnetic  resonance 
and  relaxation  for  the  ’Si  system  is  reasonably  straightforward  and 
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well  known,  partially  because  of  the  spin  l/2  nucleus  and  partially 
because  of  the  several  orders  of  magnitude  variation  in  the  Larmor  fre- 
quency, w;  the  NMR  linewidth,  = l/T^;  and  the  spin-lattice  relaxation 
rate,  l/T]. . These  quantities  are,  respectively,  ~ 5 x 10  (10,000  G), 

~ 2 x 10  (0.4  G),  and  less  than  0.01  s~ ; . The  latter  values  are  suf- 
ficiently well  removed  from  one  another  to  allow  the  choice  of  other 
laboratory  parameters , such  as  rf  amplitude  and  sweep  rates  for 
resonance  display,  which  simplify  somewhat  the  necessary  theoretical 
treatment  for  understanding  the  laboratory  measurements. 


REVIEW  OF  LITERATURE 


Since  the  pioneering  experiments  with  nuclear  magnetism  by  Rabi, 
Purcell,  and  Bloch,"  extensive  NMR  literature  has  developed,  mostly 
involving  high-resolution  spectra  of  nuclei  in  molecules  for  the  pur- 
pose of  determining  chemical  structure.  This  review  is  restricted  to 
studies  of  solids,  semiconductors,  and  in  particular  to  relaxation 
processes . 


Attempts  to  account  for  spin-lattice  relaxation  by  direct  and 
Raman  type  phonon-nuclear  interactions,  in  analogy  with  the  work  of 
Waller4  on  electron  paramagnetic  relaxation,  proved  to  be  unsuccessful 
since  the  predicted  relaxation  times,  TL , were  much  longer  than  were 
observed  experimentally.  Bloembergen5 j8  found  that  the  introduction  of 
paramagnetic  impurities  (unpaired  electrons)  produced  substantial 
changes  in  the  relaxation  rate,  1/T-l  and  its  temperature  dependence. 

He  proposed  an  order  of  magnitude  theory,  based  on  a diffusion  of  spin 
energy  to  the  lattice  via  the  paramagnetic  centers,  to  account  for  the 
increased  relaxation  rates  observed.  This  spin  diffusion  theory  has 
been  extended  by  many  authors  and  has  become  a standard  model  for  de- 
scribing nuclear  relaxation  in  solids. 


The  importance  of  paramagnetic  electrons  in  nuclear  relaxation  was 
confirmed8  through  relaxation  rate  studies  of  LiF  before  and  after 
X-ray  bombardment.  The  change  in  relaxation  rate  due  to  irradiation 
indicates  that  NMR  can  be  used  as  an  experimental  technique  to  study 
radiation  effects. 


Conduction  electrons,  interacting  with  nuclei  via  the  contact  part 
of  the  hyper fine  interaction,  can  dominate  the  relaxation  behavior  in 
semiconductors  near  room  temperature.  Theoretical  curves,  describing 
the  conduction  electron  contribution  to  the  temperature  dependence  of 
Tl  for  the  4Si  system,  were  proposed  several  years  before  any  NMR 
experimental  data  on  semiconductor  silicon  became  available. 


In  1955 » Holtzman*  and  associates  succeeded  in  measuring  Tx  at 
room  temperature  in  purified  fused  silica.  They  found  Tl  to  be  10  hours 
with  a 9-2  kG  static  field  and  1.3  minutes  in  the  Earth's  field.  Since 
the  above  initial  work,  several  authors  have  used  NMR  to  investigate  the 
J9Si  nuclear  relaxation  resulting  from  phosphorus  and  other  impurities. 


For  the  present  investigation,  phosphorus,  a substitutional  impurity 
covalently  bound  to  the  silicon  lattice,11  and  interstitial  lithium,  in 
a variety  of  complexes  with  oxygen  and  vacancies,12’13  serve  as  donors 
supplying  the  bulk  of  the  electrons  responsible  for  the  relaxation 
process . 


Room  temperature  studies14’15’16  of  phosphorus-doped  semiconducting 
silicon  have  revealed  that  nondegenerate  carriers  control  the  relaxation 
rate,  Ti-1,  for  concentrations  between  5 x 1016  and  1018  carriers/cm  \ 
For  lower  carrier  concentrations,  the  29Si  relaxation  was  attributed  to 
unidentified  paramagnetic  centers,  which  are  not  ionized  at  room  tem- 
perature. Their  effect  on  nuclear  relaxation  depends  on  neither  the 
elemental  donor  impurity  added  to  the  silicon,  nor  the  sample  prepara- 
tion. In  p-type  silicon,  Tx  was  found  to  be  approximately  eight  times 
longer  than  in  n-type  material  of  the  same  carrier  concentration 
(5  x lO'-lO18  carriers/cm3).  This  difference  becomes  smaller  with  de- 
creasing carrier  concentration.  For  mobile  carrier  concentrations  near 
1014  carriers/cm3,  the  room  temperature  relaxation  rates  in  p-type 
silicon  were  found  to  be  approximately  1.5  times  longer  than  in  n-type 
material . 


Combrisson  and  Solomon1^  measured  TL  for  the  °Si  system  in  the 
semiconducting  range  (5  x 1016  P/cm3)  at  300,  77  5 and  4.2  K but  not  at 
intermediate  temperatures.  Their  values  at  300  K (5  min)  and  77  K 
(9  min)  agree  with  the  results  of  the  present  investigation  in  a similar 
impurity  concentration  range.  They  used  dynamic  polarization  tech- 
niques18’19 to  measure  the  ?9Si  spin-lattice  relaxation  time  at  4.2  K. 


Samples  with  donor  impurity  concentrations  both  higher  and  lower 
than  5 x 1016  P/cm3  were  considered  in  the  present  investigation.  With 
decreasing  temperature  the  quantity  TL  was  found  to  increase  to  a maxi- 
mum near  50  K by  a factor  as  high  as  7-5  times  that  at  300  K and  then 
decrease  at  lower  temperatures.  In  the  same  sample,  the  equilibrium 
magnetization,  Ms,  as  measured  by  NMR  techniques  increases  with  de- 
creasing temperature  (following  Curie's  Law)  but  reaches  a maximum  near 
35  K.  At  25  K the  quantity  Ms  is  approximately  one-half  the  maximum 
value . 


As  the  donor  impurity  concentration  is  increased  above  2 x 1018 
P/cm",  the  impurity  electron  wave  functions  strongly  overlap  and  the 
electrons  become  delocalized  forming  an  impurity  band.  When  this  occurs 
the  low  temperature  electrical  resistivity  and  Hall  coefficient  measure- 
ments indicate  a transition  from  semiconductor  to  metal.20  Several 
authors  have  used  NMR  to  investigate  this  transition  in  phosphorus  doped 
silicon.  1 - ’4  Their  results  for  concentrations  less  than  1018  P/cm"  can 
be  compared  with  the  findings  of  the  present  investigation.  The  work  of 
Sundfors  and  Holcomb  is  typical  of  these  researches.21  They  measured  Tv 
for  the  9Si  system  at  1.6,  4.2,  20,  78,  and  300  K for  impurity  concen- 
trations in  the  range  from  101  to  10?c  P/cm3.  A Knight  shift  appeared 
in  the  29Si  resonance  frequency  at  the  metallic  transition.  Recent 
measurements24  at  4.2  K give  this  critical  concentration  for  impurity 
band  formation  as  L - 4 x 10 18  P/cm"  for  phosphorus-doped  samples. 


NMR  techniques  have  been  used  to  study  the  sendconductor-to-metal 
transition  in  other  solids.  Alexander25  observed  both  the  3Si  and  1 3C 
NMR  spin-lattice  relaxation  rates  in  nitrogen-doped  SiC . Look33’3  and 
associates  studied  CdO  which  was  heavily  doped  with  indium  and  observed 
nonexponential  recovery  of  the  1 1 3Cd  NMR  line  at  4.2  K.  Adams38’39  and 
associates  studied  the  semiconductor-to-metal  transition  in  chlorine- 
doped  CdS . 


In  semiconducting  silicon  at  low  temperatures,  the  nuclear  spin 
diffusion  model  can  be  used  to  define  two  classes  of  ?9Si  nuclei  with 
different  relaxation  behavior.  Nuclei  near  paramagnetic  sites  couple 
with  the  localized  electrons  via  the  scalar  contact  hyperfine  inter- 
action or  the  direct  dipole-dipole  interaction  and  proceed  quickly  to 
thermal  equilibrium  with  the  lattice.  Nuclei  farther  from  the  para- 
magnetic centers  approach  equilibrium  at  a spin  temperature  different 
from  the  lattice  temperature,  and  then  exchange  energy  with  the  lattice 
by  spin  diffusion  to  the  impurity  centers  at  a much  slower  rate . The 
latter  process  is  the  result  of  a direct  dipole-dipole  coupling  between 
more  distant  nuclei  and  the  nuclei  near  the  paramagnetic  centers . 

Jerome,  Ryter,  and  Winter3  have  presented  experimental  evidence  to 
show  that  29Si  nuclei  near  phosphorus  donor  sites  are  relaxed  by  the 
contact  interaction.  A solution  to  the  spin-diffusion  equation  for  this 
case  has  been  made  by  Buishvili.31 


The  "hopping"  model  for  impurity  conduction  of  Miller  and 
Abrahams  2 was  extended  and  applied33  to  the  relaxation  via  localized 
electrons . These  authors  considered  the  small  but  nonzero  overlapping 
of  the  electron  wave  functions  near  donor  impurities  as  an  ingredient 
of  a tunnel  phenomenon.  The  "hopping"  of  an  electron  from  donor  site 
to  ionized  donor  site  modulates  the  dipole-dipole  interaction  which  con 
trols  spin  diffusion.  A distribution  of  hopping  frequencies  implies 
that  only  those  which  satisfy  certain  criteria  will  participate  in  the 
relaxation.  Jerome  and  Winter33  were  able  to  observe  this  relaxation 
mechanism  by  controlling  the  number  of  un-ionized  donor  sites  by  compen 
sating  n-type  silicon  (6  x 10 15  P/cm3)  with  boron.  The  experiments  to 
demonstrate  this  mechanism  were  performed  at  liquid  helium  temperatures 
where  the  electron  spin-lattice  relaxation  time,  Tie,  is  large  compared 
to  a nuclear  precession. 


In  high  resistivity  silicon  (60  to  2,500  Q-cm)  TL  was  measured  as 
a function  of  field  and  temperature  using  adiabatic  fast  passage  tech- 
niques.34 TL  was  reported  to  be  independent  of  field  between  10  and 
10,000  G (290  min  ± 5%)  at  room  temperature.  In  the  temperature  range 
from  195  to  565  K the  relaxation  rate  could  be  described  by  the  empirical 
formula  Tl“1  = A exp(-W/kT),  where  W is  an  activation  energy  on  the 
order  of  0.07  eV.  An  explanation  using  the  spin  diffusion  theory  implied 
an  increase  in  the  number  of  relaxing  centers  with  increasing  temperature 
This  is  possible  if  some  of  the  impurities  are  excited  from  a diamagnetic 
state  to  a paramagnetic  state.  The  spin  system  then  relaxes  by  spin 
diffusion  to  these  paramagnetic  centers.  At  temperatures  above  565  K 
the  intrinsic  conduction  of  silicon  masks  these  effects. 


Recent  investigations16  done  in  this  laboratory  indicate  that  the 
spin  diffusion  theory  could  be  used  to  explain  an  observed  field  de- 
pendence at  room  temperature  in  moderate  resistivity  samples  (10  ft-cm). 
This  relaxation  behavior  was  attributed  to  the  formation  of  oxygen  re- 
lated paramagnetic  complexes.  In  addition,  low  temperature  data  on  a 
10  fi-cm  phosphorus  doped  sample  with  10JS  Li/cmr  added,  revealed  a two- 
fold increase  in  Tj.-1  with  decreasing  temperature  between  50  and  25  K. 

The  extensive  experimental  and  theoretical  supporting  literature 
and  the  sensitivity  to  changes  in  its  electronic  properties  make  semi- 
conducting silicon  a fruitful  system  for  investigating  radiation  damage 
effects.  Consequently,  much  effort  has  been  spent  on  its  electronic 
behavior  after  exposure  to  a variety  of  energetic  irradiations.  The 
earliest  theories  on  radiation-induced  defects  were  developed  from 
experiments35-37  on  germanium  performed  at  Purdue  University  and  at 
Oak  Ridge.  These  suggested  that  the  formation  of  vacancy- interstitial 
pairs  (Frenkel  pairs)  and  other  simple  defects  were  responsible  for  the 
large  changes  in  sample  conductivity  after  irradiation.  It  was  soon 
realized  that  the  defect  electronic  structures  are  complex  and  rely  on 
the  nature  and  concentration  of  sample  impurities . 38 s 39 

Rot  all  radiation-induced  defects  are  paramagnetic,  but  those 
which  are  can  relax  the  "Si  nuclear  system.  A comprehensive  discussion 
of  the  many  centers  which  have  been  observed  in  silicon  can  be  found  in 
the  book  by  Vook.4° 

Aside  from  their  contribution  to  nuclear  relaxation,  these  centers 
are  of  interest  for  their  role  in  semiconductor  technology.  Irradiation- 
induced  damage  sites  often  act  as  recombination  centers  which  trap 
electrons  and  shorten  the  minority  carrier  lifetime.41  This  effect  de- 
grades the  performance  of  silicon  devices  used  in  radiation  environments. 

In  1966,  J.  J.  Wysocki41  and  associates  discovered  that  the 
addition  of  lithium  to  silicon  solar  cells  improves  their  resistance  to 
radiation  damage.  They  attributed  this  effect  to  the  mobility  of 
lithium.  When  lithium  interacts  with  a damage  site,  it  may  form  a com- 
plex which  no  longer  acts  as  a recombination  center.  Thus,  lithium 
serves  to  anneal  radiation-induced  defects.  Faith4  and  associates  and 
others4 3 have  used  junction  capacitance  techniques  to  study  this  anneal- 
ing process  in  electron-irradiated  silicon  solar  cells. 


REVIEW  OF  THEORY 

Theoretical  aspects  of  nuclear  magnetic  resonance  and  relaxation 
are  reviewed  as  they  pertain  to  this  investigation.  The  classical  and 
quantum  mechanical  theories  of  weakly  interacting  spins43-45  complement 
each  other  in  illustrating  the  gross  aspects  of  NMR . Several  texts 
have  developed  these  theories  and  may  be  consulted  as  references . 4P-4 8 
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Nuclear  magnetic  resonance  of  strongly  interacting  spins,  as  it 
applies  to  the  present  investigation,  will  be  reviewed  in  terms  of  the 
recent  theory  due  to  Provotorov49  and  its  extensions  by  Goldman90’51 
and  Weber.59  Though  no  attempt  is  made  in  this  report  to  study  the  line 
shape  problem  in  silicon,  the  results  of  these  calculations  are  needed 
to  justify  the  use  of  the  adiabatic  fast  passage  technique  for  relating 
the  recorded  NMR  line  to  the  magnitude  of  the  sample  magnetization.  The 
latter  technique  has  been  used  by  most  investigators  studying  the  9Si 
system.  A notable  exception  is  the  use  of  pulse  techniques  by  Sundfors 
and  Holcomb. 9 1 


Also,  nuclear  spin  lattice  relaxation  mechanisms  will  be  reviewed. 
As  Bloembergen  and  many  others  have  demonstrated,  paramagnetic  electrons 
control  the  relaxation  rate,  TL-1,  and  the  observability  of  the  NMR 
effect  in  most  solids.  In  a semiconductor,  such  as  silicon,  impurities 
and  crystalline  defects  provide  the  paramagnetic  electron  system.  It  is 
convenient  to  divide  this  paramagnetic  concentration  into  nonlocalized 
and  localized  systems.  At  higher  temperatures  and/or  higher  donor 
impurity  concentrations,  conduction  electrons  can  produce  field- 
independent  and  field-dependent  contributions  to  nuclear  relaxation. 
These  results  from  the  interactions  of  conduction  electrons  with  99Si 
nuclei,  either  directly  with  the  P9Si  nuclei  or  indirectly  through 
electron  trapping  at  donor  impurity  sites.  The  field-dependent  contri- 
bution to  Tt-1,  due  to  the  interactions  of  conduction  electrons  with 
ionized  donor  impurity  sites,  has  received  little  attention  in  the 
experimental  literature. 


With  deeper  trap  sites  and  lower  temperatures , localized  electronic 
effects  become  important.  To  describe  their  role  in  nuclear  relaxation, 
the  interactions  between  localized  electrons  and  nuclei  are  considered. 
Also,  those  aspects  of  spin -diffusion  theory  which  are  pertinent  to  this 
investigation  are  reviewed. 


NMR  of  Weakly  Interacting  Spins 


When  a nuclear  magnetic  moment,  jjn,  interacts  with  a magnetic 
field,  H,  its  energy  is  quantized.  The  Hamiltonian  for  this  nuclear 
Zeeman  interaction  is  written 


The  quantity  pn  is  collinear  with  the  nuclear  angular  momentum, 
and  is  related  to  L by  the  gyro-magnetic  ratio  y ; 


The  simplest  quantized  system  occurs  for  a static  field,  H = H0ez. 
the  allowed  energies  are  proportional  to  the  strength  of  the  field 
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The  azimuthal  quantum  number,  m,  takes  on  (21  +1)  half-integer  values 
(m  = -I,  -I  + 1,  I - 1,  I).  For  the  9Si  nucleus  I = l/2  and 

yn  -846  Hz/G.'5”'  The  Zeeman  energies  for  J9Si  nuclei  are  sketched  on 
the  following  energy  level  diagram. 


E./; 

E-y, 


AE  = E . - E_,a  = -7nb.H0 


If  a spin  is  originally  in  the  state  |m),  a transition  to  the 
state  !m' ) will  occur  if  a perturbing  field  is  present  to  induce  the 
transition.  The  simplest  perturbing  field  is  a rotating  field  ortho- 
gonal to  the  static  field; 


Hj_  = Hi  cos  (cot)ex  + % sin  (ut)ey 


(4) 


This  field  may  be  thought  of  as  an  applied  oscillating  field  or  as  a 
spectral  component  of  the  local  field  acting  on  the  nucleus.  In  either 
case  the  interaction  is  written 

"Ki  = -y.n  • Hx  = -7n?i[HL  cos  (mt)lx  + Hx  sin  (wt)ly]  (5) 


or 


Hl  = - i y hHL[l+exp  (-iwt)  + I_exp  (icot)]  (6) 

In  Eq.  (6)  I+  and  I_  are  the  quantum  mechanical  raising  and  lowering 
operators  defined  by 


I±  = Ix  ± ily  . (7) 

First-order  time-dependent  perturbation  theory  determines  the 
probability  per  unit  time  for  transitions  between  the  states . The 
familiar  result  is  written 

Pm->m'  = (2rt/K)(7£*?Hf/4)|(m'  ll+  + I Jm)  I ?(g(m)/h  ) (8) 

where  g(w)  is  the  normalized  line-shape.  For  free  spins,  g(w)  is 
essentially  the  Dirac  delta-function  with  origin  at  the  energy  level 
separation  to  = AE/fi . The  width  of  g(to)  is  determined  by  the  transition 
time  between  the  states.  The  perturbation  expression  reveals  that 
transitions  occur  only  between  the  quantized  states  and  that  the 
resonance  response  occurs  at  the  frequency  AE/K . In  a field  of  10  kG 
the  ?9Si  NMR  occurs  at  approximately  8.46  MHz. 

Before  proceeding  to  a macroscopic  system,  it  is  instructive  to 
review  the  Heisenberg  equation  of  motion  for  a free  spin.  This 
equation  of  motion  has  the  following  form  for  an  arbitrary  field,  H: 

dL/dt  = (l/«)!Lt(]  = (7n/i*)[g  • fe,L]  . (9) 
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By  expanding  a component  of  Eq.  (9) , the  explicit  form  of  the  equation 
of  motion  is  implied; 

dLz/dt  = (7n/ih ) (Hx[Lx,Lz  ] + Hy! Ly,Lz  ]} 

/ 

= ~ln(HxLy  ‘ %-Lx) 

= 7n(L  x H)z  . (10) 

The  commutation  relations,  [L.j_,Lj]  = were  used  to  simplify 

the  terms.  By  combining  Eq.  (2)  and  (10),  the  correspondence  of  the 
quantum  mechanical  and  classical  descriptions  of  the  motion  is  demon- 
strated; 

dpn/dt  = y ndn  * 2 • ( tl ). 

This  result  is  general  but  will  be  restricted  in  further  developments 
to  a spin  l/2  system. 

In  a typical  NMR  experiment  one  deals  with  more  than  a single 
spin.  Consider  a system  of  N weakly  interacting  spin  l/2  nuclei  in  a 
static  magnetic  field.  If  this  system  were  prepared  in  a state  where 
the  individual  spins  have  equal  probability  of  occupying  either  of  their 

N 

quantized  states,  the  bulk  nuclear  magnetizatior  M = £ would  be 

i = 1 

zero.  At  normal  temperatures  and  in  a static  field,  H Hoez,  thermal 
interactions  of  these  spins  with  their  surroundings  (lattice)  result  in 
a dynamic  equilibrium  governed  by  the  Boltzmann  function, 

N./N+  = exp  (->n*Ho/kTL)  . (12) 

N±  are  the  populations  of  the  m = ± l/2  states  and  TL  is  the  temperature 
of  the  lattice. 

The  approach  to  equilibrium  is  governed  by  two  rate  processes;  i.e., 

(1)  TL ~ 1 , the  rate  of  energy  exchange  with  the  lattice,  and 

(2)  Ta-1,  the  rate  of  magnetic  energy  distribution  within 

the  spin  system. 

Experimentally  one  finds,  when  the  static  magnetic  field  is  turned 
on  or  off,  the  approach  to  equilibrium  is  governed  to  a good  approxi- 
mation by 

dMz(t)/dt  = iMo  - Mz(t)]/T1  . (13) 

Thus,  when  the  field  is  turned  on 

M*(t)  = Mo[l  - exp  (-t/TL)]  , (14) 
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and  when  the  field  is  turned  off 

Mz(t)  = Mo exp  (-t/Ti)  . (15) 

In  these  equations  Mz(t)  is  the  z-component  of  the  bulk  magnetization 
at  time  t and  is  the  equilibrium  fEq.  (l4)]  or  initial  [Eq.  (15)1 
magnetization  in  the  field  of  strength  H0 . 

The  Bloch  equations  were  written  to  describe  the  evolution  of  the 
magnetization  for  a system  of  weakly  interacting  spins . They  include 
the  torquing  effect  of  the  applied  field  and  both  relaxation  effects 
through  the  parameters  Ti_  and  T2.  Since  the  rate  process  involving  T2 
conserves  energy  within  the  spin  system,  this  relaxation  process  in- 
volves components  transverse  to  the  static  field.  Note  also  that  the 
linearity  of  Eq.  (ll)  and  (13)  is  used  in  generalizing  to  the  macro- 
scopic form; 

dMx/dt  = 7 n(M  x H)x  - (Mx/T2)  , 

dMy/dt  - 7n(M  x H)y  _ (My/T2)  > and  (Bloch  equations) 

dMg/dt  = 7n(M  x H)z  - (Mz  - MqJ/Tj.  . (16) 

In  the  presence  of  a rotating  field  in  the  form  of  Eq.  (4),  it  is  con- 
venient to  write  the  Bloch  equations  in  a reference  frame  rotating  with 
frequency  w u'ez.  This  transformation  yields 

dMx  1 / dt  7nMy  (H0  + w/^n)  " (Mx*/T2)  , 

dMy-/dt  ---  7nlMz'Hi.  “ Mx>(Hq  + w/7n)]  “ (My/T 2)  , and 

dMz'/dt  = -7nMy' H]_  - (Mz>  - Mo)/T!  , (17) 

where  Mx«,  My-,  and  Mz>  in  the  rotating  frame  correspond  to  Mx,  My, 
and  Mj  in  the  laboratory  frame . 

The  steady  state  solution  in  the  rotating  frame  (dM/dt)rot  = 0, 
is  particularly  illuminating; 

w „ m 1 (co/7  n ) + HQ]T2 

x ' M)'nHlT-  5 + .nHt'TiTa  + (l  W7nJ  + hJt2}-  ’ 


Mv-  - Mo7nH)T5  

7 ^ ~ 1 + 7 hHl  T2  + 

Mz>  = Mo  1 + f [ (^AnL. 

^ T 1 P'V  T»  T f r f . . 


7nH.  T,T,  + (;  (u/7n)  + H0  ]T a)s 


, and 


1 + {l(w/7b 

+,  H(? 

JT2r 

1 + 7 fPi'Tx  T2  + {[ 

(w/>n) 

+ H0lTo}? 

L 3 

f * 
r - 


These  equations  predict  saturation  and  line  shape  properties  for  weakly 
interacting  spin  systems;  e.g.,  in  liquids. 
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The  experimental  arrangement  detects  an  oscillating  component  of 
the  nuclear  magnetism  in  the  laboratory  frame.  For  example, 

Mx  = MX’COs(ut)  + My> sin  (wt)  . (19) 

Changes  in  ^ (or  My)  correspond  to  changes  in  the  complex  susceptibility 

x(w)  = x' (w)  + X"(w)i  (20) 


Thus , 


Mx(t)  = 2H1[x'(w)  cos  (wt)  + x"^)  sin  ( o >t ) 3 . (21) 

The  linear  driving  field  was  written  as  a rotating  field  in  Eq.  (4)  to 
simplify  the  analysis.  In  Eq.  (21)  the  amplitude  of  the  linear  field, 
2HL , is  used.  The  component  in  phase  with  the  driving  field  is  chosen 

to  be  x ' (to) ; 


'(<*>)  = (xowoTp/2) 


at. 


i + wfr^  + (at2)' 


and 


(22) 


where 

Xo  ■=  3 A = w + >nHo  j 

wo  = ?nHo  3 and  (ox  = mHi 

The  choice  of  these  susceptibility  components  follows  from  the 
effect  of  the  total  nuclear  susceptibility,  x(w)s  on  the  impedance  of 
an  LRC  tuned  circuit.  The  nuclear  system  changes  the  inductance  of  the 
circuit  from  £0  to 


(23) 


(24) 


£ = £0[l  + 4rtqx(w)  ] , (25) 

where  q is  a filling  factor.  The  susceptibility  of  the  electron  system 
is  ignored  in  Eq.  (25).  The  impedance  of  the  coil  becomes 

z = [R0  + 4irq^0wx'  (to)]  + ico/0 [ 1 + 4rtqx’(w)l  • (26) 

In  these  expressions 

X ' (to)  (Given  by  Eq.  (22)]  reflects  changes  in  the  frequency 

of  the  tuned  circuit  with  no  net  absorption  or  emission 
of  rf  energy.  The  component  x'(0  is  called  the  dis- 
persion mode  of  the  NMR  signal. 


— • ^ 


X"(w)  [Given  by  Eq.  (23)]  reflects  a net  absorption  or 

emission  of  rf  energy  by  the  sample.  The  component 
X”(w)  is  called  the  absorption  mode. 

Phase  detection  techniques  may  be  used  to  select  either  of  these  modes 
or  a mixture  of  these  modes. 

NMR  in  Solids 

The  Bloch  equations  adequately  describe  the  macroscopic  properties 
of  NMR  for  gases  and  many  liquids,  where  Tt  = T2.  For  strongly  coupled 
spins  as  in  a solid,  To  « Tx . In  these  solids,  when  the  rf  level  is 
small  compared  to  the  spin-lattice  interaction  (ynHfliTo  « l),  the 
populations  of  the  nuclear  Zeeman  levels  remain  unchanged  during  the 
measurement  and  NMR  lineshapes  similar  to  those  predicted  by  the  Bloch 
equations  are  observed.  However,  for  larger  rf  levels  (7nHiTi.T2  » l) 
the  coupling  of  the  nuclear  Zeeman  system  with  the  rf  field  must  be 
included.  Experimental  results54  on  the  absorption  and  dispersion  modes 
of  copper  and  aluminum  indicate  that  the  dispersion  signal  does  not 
saturate  in  the  manner  predicted  by  the  Bloch  equations  and  saturational 
narrowing  is  observed.  Redfield54  used  thermodynamic  arguments  to 
modify  the  Bloch  equations  for  rf  levels  which  are  large  compared  to  the 
low  field  dipolar  line  width  (Hi.  » D/Vn)  • This  theory  has  been  de- 
veloped and  confirmed  by  experiment.55-58 

Goldman59**  and  associates  have  used  this  thermodynamic  approach 
to  describe  the  NMR  dispersion  derivative  signal,  which  is  recorded 
during  an  adiabatic  fast  passage.  Such  a passage  is  defined  by  the 
inequalities 

1/Tj.  « A/D  « 1/Ta  , (27) 

where  A Vn  is  the  sweep  rate  (G/s)  for  the  static  field  and  D/yn  is  the 
steady-state  linewidth  in  the  low  rf  field  limit . In  '^he  present  in- 
vestigation the  Inequalities  (27)  are  satisfied;  Ti  > 100  s,  To  < 1 ms, 
and  a/D  ~ 1 s- • . The  components  Mxi  and  Mz«  in  the  rotating  frame  are 
given  by 


Mx<  oc 


Ta- 


+ 


~ • n \ r 

D-r 


and 


Mz, 


A 

(t'f"  + wy> 


(28) 


Equations  (28)  imply  that  the  magnetization  remains  parallel  to  the 
effective  field  during  the  adiabatic  fast  passage;  components  of  the 
magnetization  orthogonal  to  the  effective  field  decay  to  zero  in  a time 


*Goldman,  Ref.  51. • pp.  41-43. 
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of  the  order  to  T2,  and  after  completion  of  the  passage  the  total  mag 
netization  is  reversed.  The  dispersion  derivative  signal,  which  is 
recorded  in  the  experiment,  is  observed  to  be  asymmetrical  because  of 
partial  saturation  is  not  predicted  by  Eq.  (28)  and  is  left  to  a 
solution  of  the  Provotorov  equations  and  their  generalizations. 


The  Provotorov  Equations 


In  the  frame  of  reference  rotating  with  the  applied  rf  field 
Hamiltonian  for  the  spin  system  is  written 


Here,  H is  the  secular  part  ([Hs,  Iz]  = 0)  of  the  nuclear  spin-spin 
interaction  in  the  rotating  frame,  A r m - w0,  and  Ui  = 7nHi • Oscillat- 
ing components  and  interactions  with  the  lattice  are  not  included  in 
Eq.  (29).  For  small  rf  fields  (Hl  « Hd)  and  far  from  the  center  of  the 
resonance  line  (a  » nHd),  the  Zeeman  and  dipolar  terms  may  be  con- 
sidered separately  as  constants  of  the  motion.  The  second  condition 
may  be  relaxed  through  a development  of  the  Provotorov  equations  in  the 
laboratory  frame.*  Provotorov  suggested  that  the  Zeeman  and  dipolar 
systems  come  to  equilibrium  separately  in  a time  on  the  order  to  T2  and 
can  thus  be  described  by  separate  reciprocal  temperatures  (01,0).  The 
density  matrix  formalism  is  used  to  describe  the  system  for  T2  « t « T 
The  density  matrix  is  written 


where 


and  <72  contains  the  coupling  term  f which  produces  an  evolution  of 
a and  0. 


In  the  interaction  representation  a second-order  perturbation 
calculation  gives  the  evolution  of  the  density  matrix  and  the  Provotorov 
equations  for  (a,0)**; 


In  these  equations  W W(a)  itos  g(A)  and  the  quantity  g(A)  is  the 
normalized  line  shape  for  the  absorption  signal  at  low  rf  levels. 


♦Goldman,  Ref.  51,  pp.  102-107 
♦♦Goldman,  Ref.  51,  pp.  75-83- 


/ dAg(A)  1 


Alternatively,  the  quantity  g(A)  is  defined  as  the  cosine  Fourier 
transform  of  the  free-induction  decay  signal.  The  relaxation  terms, 

(a  - 0!i)/T1z  and  (3  - 32)/Tp,  are  added  formally  to  allow  for  spin- 
lattice  relaxation.  Further  considerations  described  by  Goldman  give 

MX'  = -rnTr(lz)fwi3  + «u>ig'  (a)a(a  - 3)1  , 

My-  = -(W/uj.  )7nATr(lz)(a  - 3)  , and 

Mz,  - -aAynTr(lz)  . (32) 

In  these  equations  Mx> , My< , and  Mz-  are  the  rectilinear  components  of 
the  magnetization  in  the  rotating  frame.  The  quantity  g' (a)  is  the 
sine  Fourier  transform  of  the  free-induction  decay  signal.  To  complete 
the  system  of  equations,  a quantity  with  the  dimensions  of  a magnetiza- 
tion is  defined  by 


where 


Mp  = -3D>nTr(l|)  -7n((Hs)/D) 


XrMJ 

. Tr(Mz) 


The  Provotorov  equations  for  the  magnetization  in  the  rotating  frame 
thus  become 

dMz-/dt  = -W[MZ  - (a/D )Mp ] - (Mz-  - Mq)/T1z  , 

dMp/dt  W (a/D) [Mg.  - (A/D)Mp]  - (Mp  - Mp  )/Td  , 

Mx ' = (wi/DjMp  + itWig'(A)fMz-  - (a/D )Mp  ] , and 

My,  = (W/o>1)[Mz.  - (a/d)Md]  . (3M 

V/eber's  Modification  of  the 
Provotorov  Equations5 " 

In  Weber's  development  of  the  Provotorov  equations,  the  quantity 
H is  separated  into  two  commuting  parts. 


P i 


Hs  - Hss 


'The  quantities  Hss  and  Hsz  are  the  scalar  (isotropic)  components  and  the 
remaining  bilinear  z components  of  the  dipolar  field,  respectively. 

Weber  has  demonstrated  that  Hsz  is  responsible  for  the  free  induction 
decay  and  appears  in  the  density  matrix  describing  the  evolution  of  the 
nuclear  magnetization.  His  modified  Provotorov  equations  were  obtained 


r 


\ 


i 


V. '} 

i;  I 

4 i 

L 

V-  3 

h 

n 


without  the  introduction  of  the  spin  temperatures  (a, (3),  and  are  valid 
for  the  low,  intermediate,  and  high  rf  level  cases.  These  equations 
are  written  for  moderate  sweep  rates  as 

dMZ’/<it  = -W[MZ-  - (a/wd)md1  ~ (MZ’  - M0)/Tlz  , 

dMp/dt  = [top/(u£  + wf))  ](A/wp)fMz'  - (A/wp)Mp]  - (Mp  - Mp0)/Tp  > 


Mv 


= (wi/wd)Md  + Jt^g'CAjtMz'  - (A/wp)Mp]  , and 


My.  = (W/wiJtMz-  - (a/wb)Md]  . 

Here,  the  quantities  M^  and  are  similar  to  the  corresponding 
quantities  Mp  and  D in  the  Provotorov  theory.  They  are  defined  by 


(36) 


Md  = ~y  n(Hsz  )/tJD 


and 


(37) 


The  modifications  to  the  Provotorov  equations,  which  resulted  from 
Weber's  development,  are  seen  by  comparing  Eq.  (33)  and  (34)  with 
Eq.  (37)  and  (36).  Weber1 1 ’’  has  numerically  integrated  Eq.  (36)  for 
moderate  sweep  rates  and  a broad  range  of  rf  levels.  The  line  shapes 
obtained  from  these  integrations,  corresponding  to  the  dispersion  de- 
rivative signals,  exhibit  partial  saturation  and  agree  with  the  experi- 
mental line  shapes  observed  in  the  present  investigation. 

Relaxation  Mechanisms  in  Silicon 

The  relaxation  of  ~3Si  nuclei  in  a static  field,  H0et,  can  be 
effected  by  numerous  mechanisms.  Inherent  in  many  of  the  models,  de- 
scribing the  relaxation,  is  a time- dependent  local  field,  hj^t).  The 
power  spectrum  of  this  field  contributes  to  the  NMR  relaxation  rates 
Tl-1  and  T2-1,  and  shifts  the  nuclear  Larmor  frequency. 


t>y 


The  correlation  function  for  a component  of  this  field  is  defined 


f(t)  = (hL*(t)hL(t  + t)>  , 


(38) 


where  the  average  is  taken  over  all  time  for  a particular  nucleus.  The 
spectrum,  j(to),  is  obtained  through  a Fourier  expansion  of  the  correla- 
tion function; 


f ( t ) = / j(w)  exp  (icoi)dw 
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and 


J(w)  l/2« 


f f(i)  exp 
-'-00 


( -iun )di 


(39) 


The  saturation  of  NMR  lines  with  rf  fields  and  time-dependent  perturba- 
tion theory  indicate  that  spectral  components  of  h^(t)  near  the  nuclear 
Larmor  frequency  will  be  responsible  for  nuclear  relaxation.  Thus,  by 
using  the  two  equal  random  field  components  orthogonal  to  the  static 
field,  H0ez,  the  resulting  expression  for  the  relaxation  rate  at  a 
position,  r,  from  a paramagnetic  center  is 

1/Ti(r)  = 2rgj£(u>)  . (40) 

Because  of  the  dependence  on  r in  this  expression,  further  considerations 
are  needed  to  obtain  the  relaxation  rate,  T^-1,  for  the  bulk  sample. 
Several  mechanisms  affecting  t^ie  quantity  Tj.-1 , will  be  considered  in 
the  remaining  sections.  These  models  do  not  describe  the  only  possible 
relaxation  mechanisms  in  solids,  but  represent  those  which  are  important 
to  the  present  investigation. 

At  the  basis  of  all  these  mechanisms  is  the  hyperfine  interaction 
between  the  magnetic  moments  of  an  electron,  qe,  and  a 29Si  nucleus, 

Mjj.  The  interaction  is  written  in  operator  form  as 

‘Khf  = (-8*/3)ue  • tin6(r)  + [)ie  • tin  - 3(jie  • *)(|in  • ^)]/r3  . (4l) 

In  this  expression  r rr  is  the  position  vector  locating  the  nucleus 
with  respect  to  the  electron.  The  first  term  is  called  the  Fermi  con- 
tact term,  and  the  second  term  the  direct  dipole-dipole  interaction. 

Since  the  interaction  energy  can  be  written  as  a field  perturbing  a 
nuclear  dipole  £e  “[^n  ‘ matrix  elements  for  the  local  field  can 

be  written 

(hL)hf  = (ar/3)re^f(o)sn(o) 

-f  Jff(r){(ye«/r3)[S  - 3r (S  • r ) ] }\Ki(r )d3r  , (42) 

where  the  magnetic  moment  of  the  electron  has  been  written  in  terms  of 
its  angular  momentum,  = y ehS . The  quantities  \|fj_(r)  and  >Kf(r)  are 
the  initial  and  final  electronic  states,  respectively. 

Often  the  initial  and  final  electron  states  are  the  same  spatially. 
In  this  case  the  matrix  elements  [Eq.  (42)]  represent  the  local  field 
and  the  contact  term  is  proportional  to  the  quantity  |l|/(0)|  • The  non- 
localized  electron  density  is  written  l^(0)l  = va  ) ( 0 ) I , where 

^k(o)  is  averaged  over  the  Fermi  surface  and  is  normalized  to  the  atomic 
volume.  For  localized  electrons  the  ground-state  wave  function  is  also 
influenced  by  the  periodicity  of  the  silicon  lattice.  Descriptions  of 
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the  shallow-donor  ground-state  wave functions  in  silicon  have  been 
proposed  by  Kohn  and  Luttinger61 »6t  and  recently  by  Ivey  and  Mieher.  3 

The  degree  of  electron  localization  can  have  a marked  effect  on 
the  relaxation  rate  for  the  ‘?9Si  system.  It  is  thus  natural  to  divide 
a discussion  of  relaxation  mechanisms  in  semiconducting  silicon  into 
nonlocalized  and  localized  electron  effects.  The  transition  region  will 
also  be  considered. 

Nuclear  Relaxation  Due  to 

Conduction  Electrons 

The  statistical  behavior  of  the  conduction  electrons  largely 
controls  the  temperature  dependence  of  the  nuclear  relaxation  rate. 
Korringa64  has  developed  a general  expression  for  the  contribution  to 
nuclear  relaxation  due  to  conduction  electrons . This  relaxation  rate  > 
is 


1/TX)C  = (2/Vi ) [v^/  (2tt  )5  ] f dk  f dk'  |e(kk* ) |?[l  - F(E+)] 

x F(E+  + guBH0/kT)&(Ev  - e£  - gUBH0)  , (43) 

where  the  unprimed  and  primed  quantities  represent  initial  and  final 
electron  states,  respectively.  In  this  expression,  the  energy  of  the 
electron  state,  (k,S  = ± l/2),  is  E£  and  the  interaction  which  effects 
nuclear  transitions  is  the  Fermi  contact  interaction; 

= (S+I.  + S_I+)e(kk'  ) 

and 

e(kk’)=  (-8n/3)tf:37e>ncp£- (0)^(0)  . (44) 

In  silicon,  with  conduction  electron  concentrations  less  than 
roughtly  1017  electrons/cm3,  the  electron  system  is  nondegenerate615 
and  the  Fermi  function  F(|)  reduces  to  the  Maxwell-Boltzmann  distribu- 
tion. Bloembergen  has  applied  Eq.  (43)  to  the  nondegenerate  electron 
system  and  has  obtained  the  following  expression  for  the  relaxation  rate 

1/Tl)c  = (512/9)vfrl7nNo(|cPk(0)|;?)2(2nm!|m|kT)'/j  . (45) 

The  quantities  N0,  mu,  and  m^  are  the  nondegenerate  conduction  electron 
density  and  the  longitudinal  and  transverse  effective  mass  components, 
respectively.  Expression  (45)  has  been  modified  to  include  the  six  con- 
duction band  minima  in  the  silicon  lattice  and  was  fitted  to  room  tem- 
perature experimental  data  to  obtain  vf  ( ('Dk(O)  \ 2Y '. 14  ,ls 

It  has  been  tacitly  assumed  in  the  previous  applications  to  semi- 
conductors, that  those  electrons  not  contributing  to  nuclear  relaxation 
by  direct  electron-nuclear  interactions  are  ineffective  in  relaxing  the 
nuclear  system.  However,  through  the  quantity  |e(kk')|  in  Eq.  (43), 


modulations  of  the  local  field  produced  by  electron  trapping  at  donor 
sites  can  be  included  as  a relaxation  mechanism.  The  effects  of 
trapping  will  appear  as  a frequency  (field)  dependence  in  the  nuclear 
relaxation  rate,  if  the  power  spectrum  of  the  local  field  is  nonuniform 
in  the  range  of  the  nuclear  Larmor  frequency. 

Middleton'6  describes  a wealth  of  analytical  techniques,  which 
have  been  applied  to  such  stochastic  processes . The  power  spectra  re- 
sulting from  many  of  these  calculations  contain  a frequency  dependence 
in  the  range  of  the  NMR  frequencies  used  in  the  present  investigation. 
Two  simple  spectra  include  the  collision  frequency  and  residence  time 
for  an  electron  at  a donor  site. 

Suppose  the  local  field  at  a nucleus  consists  of  a series  of 
rectangular  field  pulses  due  to  trapping  of  electrons  at  nearby  para- 
magnetic sites.  The  pulses  are  assumed  to  appear  with  average  frequency 
A and  have  an  average  duration  t^.  If  the  occurrence  of  these  pulses 
obeys  a Poisson  distribution  and  a distribution  of  pulse  heights  is 
allowed,  the  power  spectrum  for  the  local  field  has  the  form 

JM  = OC,b)'»W  * ^ • <w> 

The  quantity  Tb  is  proportional  to  the  residence  time,  but  is  shorter. 

It  represents  the  time  the  local  field  assumes  a particular  magnitude. 

Van  Vliet  and  Fassett6  considered  the  frequency  spectrum  for  a 
simple  generation  recombination  (g-r)  process,  and  found  the  magnitude 
of  the  local  field  spectrum  proportional  to  the  mean  square  of  the 
fluctuations  in  the  carrier  concentration; 

J(u)  «W2)/2*)[  /(I  + w?t?)]  , (47) 

where  t is  the  carrier  lifetime.  In  terms  of  the  NMR  experiment  t'1  is 
the  mean  frequency  for  modulations  in  the  electron  density  at  a nuclear 
site.  Van  Vliet  and  Fassett  extend  their  trapping  theory  to  the  shot 
noise  problem  and  to  the  stochastic  processes  in  multilevel  semiconduc- 
tors. The  trapping  spectrum  has  been  obtained  in  germanium  with 
group  III  and  V impurities  in  the  liquid  helium  temperature  range.  ' 

Relaxation  by  Localized  Electrons 

In  semiconducting  silicon  at  low  temperatures,  the  presence  of 
unpaired  localized  electrons  can  dominate  the  nuclear  relaxation  due 
to  conduction  electrons . Relaxation  by  localized  electrons  in  solids 
is  often  described  in  terms  of  the  spin  diffusion  theory  originally 
proposed  by  Bloembergen »"  In  this  model,  the  nuclear  system  is  con- 
veniently divided  into  two  parts;  i.e.,  (l)  nuclei  near  paramagnetic 
sites  are  directly  relaxed  by  the  localized  electronic  moments  through 
the  h.yperfine  interaction,  and  (2)  nuclei  farther  from  the  paramagnetic 
sites  are  relaxed  through  nuclear  spin-spin  exchange  (spin -diffusion) 


f ' 


, - 1 


to  these  sites. 


The  development  and  experimental  confirmation  of  this 
theory  is  well  represented  in  the  reviews  by  Khutsishvili • > Here, 

aspects  of  the  spin-diffusion  theory  will  be  included  which  apply  to 
moderate  resistivity  silicon. 


For  a very  small  paramagnetic  concentration,  each  9Si  nucleus 
will  be  strongly  influenced,  at  most,  by  a single  paramagnetic  electron. 
The  spectrum  of  the  local  field  can  generally  be  described  by  the  cor- 
relation time,  tc.  For  the  case  considered  here  tc  will  be  determined 
by  the  electron  spin-span  relaxation  time  and/or  the  electron  spin- 
lattice  relaxation  time. 


f'l 

i ’r  v 


The  correlation  time,  tc,  and  the  rms  amplitude  of  the  local  field, 
( I hp,(t) ! ?}'/s,  determine  the  nuclear  relaxation  rate  and  the  observability 
of  the  nuclear  resonance.’1  If  7n^|hL(t)|  )'/’tc  » 1,  the  “’9Si  nucleus 
will  interact  with  the  time -dependent  magnetic  field  of  the  electron 
moment  ue ; 


(lhL(t)|?)v»  - *>e[ S(S  + l)/3]v,(l/r  ) - *7e/2r 


(48) 


In  this  approximation,  a nucleus  25  A from  a paramagnetic  center  will 
be  influenced  by  a 0.6  G field  if  xc  > 0.16  ms.  The  contribution  of 
this  nuclear  spin  would  be  lost  from  the  NMR  line.  On  the  other  hand, 
if  n ^ | hi ( "t ) | )'/itc  « 1,  the  nucleus  will  be  influenced  by  the  much 
weaker  time  average  value  of  the  electron  moment  in  the  static  field 
Hoez ; 


(heZ)  Mez  (>JezH-o/kT ) 


(49) 


where 


'ez 


b>e[S(S  + l)/3]'A  . 


Rorschach7 ~ calculated  the  spectral  density  for  the  z- component  of 
the  unpaired,  localized  electron's  magnetic  moment.  An  exponentially 
decaying  correlation  function,  with  residual  statistical  correlation 
for  times  long  compared  to  tc,  was  assumed.  His  results  are  applicable 
to  the  case  7n(  Ihp^t)  I ?)'/:,tc  « 1.  The  spectral  density  for  the  moment 
then  has  the  form 


uez' 


(w)  - (uez)?&(w)  + ((u|z)  - (nez)?)[Tc/(l  + U)  t|)]/it  , (50) 


or  explicitly  for  the  case  yehS (Ho/kT)  « 1,, 


xez 


(w)  = y4r|([S(S  + l)/3][H0/kT])  &(w) 

+ S(S  + 1)/3][tc/(1  + ^t6)1A  , 


(51) 


Since  - uex  + ney  + UgZ,  the  spectral  density,  Ju  (w),  can  be 

used  to  describe  the  low-temperature  relaxation  of  the  ?qSi  nuclei  near 
the  paramagnetic  sites; 
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If  the  direct  dipole-dipole  interaction  is  assumed  and  an  average  over 
angles  is  made,  the  relaxation  rate  near  the  sites  is  given  by 


1/Ti(r)  = c/r6  , (5 

where 

C = (2/5)*27!7nS(S  + 1)[tc/(1  + w2t2) ] . 

In  analogy  with  Eq.  (53)  the  contribution  to  l/T((r),  determined 
by  the  contact  interaction,  can  be  written 


l/Ti'Cc)  = (l6It/27^27l7^[S( S + l)/3lU(r)|4('Vc/(1  + W?T§) ] . (54) 


Buishvili  1 used  the  hydrogen-like  s-state  approximation  to  describe 
l/Ti(r)  in  his  solution  of  the  spin  diffusion  problem. 

[l/li'(r)  = const,  x exp  (a/r),  where  a is  the  Bohr  radius  of  the 
localized  electron's  wave  function.] 


Spectral  components  of  the  local  field,  approximately  in  the  range 
2jr/Ta  < u < 2rc/To,  can  "quench"  (i.e.,  remove  nuclear  spin  contributions 
from)  the  NMR  line/2  The  quenching  field  calculated  from  Eq.  (51)  is 
given  by 


/ dmg(w)JWez(w) 

127l[S(S  + l)/3]j« 


quench 


where  g(u)  is  the  nuclear  dipolar  broadened,  low-field  absorption  line 
shape.  However,  if  xc  » 7n1(|h]j(t)  |2)“'/j,  the  quenching  field  is 
determined  by  Eq.  (48). 7 1 


For  localized  electrons  tc  is  determined  by  the  smaller  of  the 
electron  spin-lattice  relaxation  time,  Tiej  and  the  spin-spin  relaxation 
time,  T2e.  The  quantity  Tle  depends  on  temperature  and  field,  but  not 
strongly  on  the  local  concentration  of  paramagnetic  centers . On  the 
contrary,  Toe,  does  not  strongly  depend  on  temperature  and  field,  but 
is  sensitive  to  the  local  concentration  of  paramagnetic  impurities. 


In  moderate  resistivity  silicon,  the  formation  of  pairs  or  clusters 
of  impurities  can  affect  the  nuclear  relaxation  rate.  Below  4.2  K and 
with  donor  impurity  concentrations  less  than  5 x 101T  P/cm',  Jerome 
noticed  an  increase  in  T}.'1  with  increasing  impurity  concentration.  He 
attributed  this  increase  in  relaxation  rate  to  the  tunneling  of  elec- 
trons from  nonionized  to  ionized  donor  sites.  This  mechanism  required 
the  silicon  sample  to  be  partially  compensated  (contain  both  acceptor 


and  donor  impurities).  The  impurity  conduction  model2 3  » > 74  was  used 
to  obtain  a contribution  to  the  nuclear  relaxation  rate  for  nuclei  near 
a pair  of  donor  sites.’3  This  relaxation  rate  is  written 


l/T{(r)  = (c/r'Mexp  (-A/kT)/[l  + exp  (-A/kT)  ];  } , (56) 

where  C is  given  by  Eq.  (53)  and  1/tc  is  interpreted  for  this  mechanism 
as  the  "hopping"  frequency.  The  quantity  A is  the  electrostatic  energy 
difference  between  the  donor  sites  due  to  the  nearest  ionized  acceptor 
center.  Approximately  5$  compensation  is  expected  in  moderate  resis- 
tivity Czochralski  grown  n-type  silicon.30 

Buishvili  8>  " and  associates  considered  the  effect  of  electron 
exchange  interactions  on  nuclear  relaxation  in  terms  of  the  spin  tem- 
perature model.  They  define  an  "exchange  reservoir"  with  inverse 
temperature  and  develop  relaxation  equations  involving  this  system. 
An  application  to  F-center  pairs  and  clusters  in  KC1  was  suggested  in 
their  references.  'r 


The  excitation  of  donor  pairs  from  the  diamagnetic  ground  state  to 
the  paramagnetic  triplet  states  provide  spectral  components  to  the 
local  field.  The  formation  of  donor  pairs  has  been  considered  by 
several  authors.  8-80  For  the  Li-doped  samples  used  in  the  present  in- 
vestigation, Jayapandian81  has  observed  a loss  in  the  integrated  para- 
magnetism of  the  donor  system  at  low  temperatures  and  he  has  attributed 
this  effect  to  the  formation  of  donor  pairs  (or  clusters). 


All  of  these  models  c.-n  provide  spectral  components  to  the  local 
field  influencing  9Si  nuclei.  Since  the  interaction  between  para- 
magnetic electrons  and  nuclei  has  a strong  radial  dependence,  the 
electron-nuclear  interaction  may  become  small  compared  to  the  nuclear 
dipolar  interaction.  In  these  low  paramagnetic  concentration  solids, 
a spin  diffusion  mechanism  is  introduced  to  account  for  the  single 
relaxation  rate  observed. 


Review  of  Spin  Diffusion  Theories 


Local  magnetization,  created  near  the  paramagnetic  sites  (or 
clusters)  by  energy  exchange  with  the  lattice,  is  transported  to  dis- 
tant nuclei  by  nuclear  spin-spin  exchange.  Field  gradients  near  the 
paramagnetic  centers  provide  a diffusion  barrier  to  spin- spin  exchange 
Three  cases  of  spin  diffusion  theory  have  received  considerable  atten- 
tion in  the  literature.  These  are 


(l)  the  "diffusion  limited"  case,  which  is  limited  by  spin 
spin  exchange; 


(2)  the  "rapid  diffusion"  case,  which  is  limited  by  the 

diffusion  barrier;  and 
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(3)  the  "diffusion-vanishing"  case,  where  several  impurities 
strongly  affect  each  nucleus . 

Cases  (l)  and  (2)  will  be  reviewed  first.  If  each  9Si  nucleus  is 
assumed  to  be  influenced  by  a single  localized  electronic  center,  the 
transport  equation  for  the  magnetization  density,  m(r,t),  including  the 
paramagnetic  sources  becomes 

W?m(r,t)  - (d/dt)ra(r,t)  = [l/Ti  (r)  ]fm(r,t)  - md  . (57) 

Here,  D controls  the  diffusion  rate  (D  a Ta-1),  l/Tx(r)  describes  the 
relaxation  behavior  produced  by  the  paramagnetic  sites,  and  niQ 
represents  the  equilibrium  magnetization  density  in  the  sample.  As 
stated  above,  l/T[(r)  will  always  contain  a substantial  radial  depend- 
ence. This  radial  dependence  must  be  included  in  a solution  of  Eq.  (57). 
The  diffusion  equation  has  been  solved  when  l/T{ (r)  is  determined  by  the 
contact  interaction  1 and  the  direct  dipole-dipole  interaction.  The 
solutions  have  similar  properties  and  can  be  discussed  in  terms  of 
several  radii  measured  from  the  sites.  These  are  listed  below  with 
approximate  expressions  for  their  size. 

1.  a is  the  Bohr  radius  in  silicon,  a - a0k(m/m*)>  where  a0 
is  the  free  hydrogen  Bohr  radius,  k is  the  dielectric  con- 
stant of  silicon,  and  m*  is  the  effective  mass  of  the 
electron . 

2.  b is  the  diffusion  barrier  radius.  At  roughly  this  radius 
the  electron's  field  differs  between  adjacent  nuclear 
sites  by  the  linewidth  of  the  nuclear  dipolar  system.  It 
is  a barrier  to  spin-spin  exchange  processes.  An  estimate 
of  this  radius  is  given  by 

b = a/2  to  [(l6/3)((uez)/Mn)]  , (58) 

where  the  contact  interaction  dominates,  and 

b - a[3(uez)/nn)1/4  , (59) 

where  the  direct  dipole-dipole  interaction  dominates. 


3.  i is  the  linewidth  barrier  radius.  Inside  this  radius 
the  field  of  the  localized  electronic  moment  is  larger 
than  the  dipolar  linewidth  for  the  nuclear  system.  An 
estimate  of  l is  given  by 

t " a/2  ^ (8/3)  ( (uez )/tJn) ) > 

where  the  contact  interaction  dominates  and 

a[(qe7)/un]1  , 
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where  the  direct  dipole-dipole  interaction  dominates. 

4.  3 = (c/d)1/4  is  the  pseudo-potential  radius  or  scattering 
length.  An  estimate  of  8 is  easily  obtained  when  the 
direct  dipole-dipole  interaction  dominates.*  Without 
spin-diffusion,  a nuclear  moment  a distance  r from  a para- 
magnetic impurity  will  be  relaxed  in  a time  tj.  r'  /C  • 

If  diffusion  is  allowed  and  the  impurity  is  removed,  a 
spin  at  a distance  r from  the  impurity  will  flip  in  a 
time  t2  = r?/D.  These  times  are  equal  when 
r » (c/D)1/4  =3.  C is  a measure  of  the  strength  of  the 
field  produced  by  the  impurity  and  depends  on  temperature 
and  possibly  the  external  field.  The  quantity  D is  the 
diffusion  constant  and  is  determined  by  nuclear  spin- 
spin  exchange  processes. 

5.  R is  the  mean  radius  of  the  volume  occupied  by  a para- 
magnetic impurity; 


R = (3/4nN)1/3  , (62) 

where  N is  the  impurity  concentration. 

The  spin  diffusion  problem  was  first  solved  analytically  by 
Khutsishvili32  and  DeGennes83  for  the  diffusion  limited  case. 

Blumberg84  considered  the  rapid  diffusion  case  and  found  non exponential 
recovery  (Mz  ~ t1/2)  for  a short  time  following  saturation.  The  more 
general  diffusion  problem  was  solved  by  Rorschach  3 for  the  direct 
dipole-dipole  interaction  and  gives  the  nuclear  relaxation  rate  for  the 
sample , 


l/Tj.  8irND3[r  (3/4 )/r  (i/4)  ][ I-3/4 (&)/i3/4(s)  ] , (63) 

where  Im(x)  i“mJm(ix)  is  the  modified  Bessel  function,  r(x)  is  the 
gamma  function,  and  5 = 32 /2b'? . 

The  limiting  behavior  of  Eq.  (63)  is  of  interest; 

rapid-diffusion  case  (5  « l)  l/Ti  -»  (8jt/3)NI)3(3  /2b  ' ) (64) 

and 

diffusion-limited  case  (8  » l)  l/Ti  -»  (8it/3)ND3  (65) 

Equations  (64)  and  (65)  lead  to  the  following  contributions  to  Tj_-1  if 
tc  « >n1<|hL(t)|2>"1/?  and  StfyeHo  « kT: 

(5  « 1)  l/Ti  a (N/b3)(Tc/fl  + (TcynH0)?]l  (66) 

and 

(6  » 1)  1/T,.  a N{tc/[1  + (tcrnHo)2])1/4  . (67) 


♦Goldman,  Ref.  51,  PP-  65-66. 


These  expressions  are  field  (frequency)  and  possibly  temperature 
dependent  (through  tc). 

Fukushima  and  Uehling8,3  solved  the  diffusion  equation  for  larger 
impurity  concentrations,  where  several  paramagnetic  centers  influence 
each  nucleus.  They  established  a criterion  for  the  diffusion-limited 
case  (p/R  « 0.37) , and  write  the  following  expression  for  larger 
impurity  concentrations : 

1/Tl  a n^pV^PdU-p)  . (68) 

In  this  expression,  p = 0.68((3/r)  + 1/4. 

In  his  thesis,  Tse36  considered  the  case  of  larger  impurity  con- 
centrations with  5 » 1 and  b < R.  He  finds 

1/Tl  (4rtND0:?/R)[K3/4(A)/KI/4(A)]  , (69) 

P /2R  and  Km(x)  (n/2)im+1Hm^ 1 ^ (ix)  is  the  modified  Hankel 
In  the  low  concentration  limit  this  expression  reduces  to  the 
limited  result  and  with  higher  concentrations  the  diffusion 
case  is  obtained; 

l/Tt  = 4rtND02/R  . (70) 

Lowe  and  Tse  in  a later  paper87  considered  the  diffusion  vanishing 
case  using  a different  approach.  Recognizing  that  l/TL'(r)  has  minima 
in  the  solid  far  from  the  paramagnetic  sites,  they  expanded  l/T{(r) 
about  a typical  minimum.  This  expansion  to  the  harmonic  approximation 
allowed  the  diffusion  equation  to  be  separated  and  the  time-dependent 
magnetization  to  be  obtained. 

In  recent  years  the  interaction  between  electron  spins  and  the 
effect  of  this  interaction  on  the  form  of  the  correlation  function  have 
become  of  interest.  This  interest  has  been  motivated  by  the  theoretical 
work  of  Blume  and  Hubbard88  and  the  NMR  spin-lattice  relaxation  studies 
of  the  intermetallic  compound  La;L_cGdcAl2 ,89  The  correlation  function 
for  the  electronic  moment  was  found  to  deviate  from  the  simple  decaying 
exponential  used  to  obtain  Eq.  (50) .3°  However,  with  the  low  donor 
impurity  concentrations  considered  in  this  experiment,  and  the  com- 
plexity of  the  impurity  system,  these  effects  can  probably  not  be 
resolved. 


where  A = 
function, 
diffusion 
vanishing 


EXPERIMENTAL  PROCEDURE 

The  NMR  spectrometer  was  a fixed-frequency,  crossed-coil  type, 
generally  similar  to  commercial  models  but  with  some  modifications. 

The  rf  transmitter  output  to  the  probe  was  capable  of  higher  power  (to 
100  W)  and  was  electronically  stabilized  in  amplitude  to  ensure  constant 
rf  field  at  the  sample  over  long  periods  of  time.  The  receiver  employed 
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rf  phase  detection  in  order  to  stabilize  selection  of  the  dispersion 
mode  and  to  increase  signal-to-noise  ratio.  The  rf  signals  were  de- 
rived from  crystal-controlled  oscillators. 


The  external  magnetic  field  was  supplied  by  an  electromagnet 
characterized  by  stability  and  homogeneity  at  the  sample  of  less  than 
25  mG  in  10  kG.  The  field  was  modulated  at  about  40  Hz  to  produce  an 
af  signal  output  from  the  receiver  which  was  further  amplified  and 
phase-detected  at  40  Hz.  The  resulting  analog  signal,  the  derivative 
of  the  dispersion  line,  was  processed  by  an  analog-to-digital  converter 
and  recorded  in  digital  form  in  successive  addresses  of  the  Nicolet 
Averager  as  the  external  field  was  scanned  over  a 50  G range  centered 
at  the  Larmor  resonance  condition.  The  recorded  data  were  then  processed 
by  the  Digital  Equipment  Corp.  minicomputer  to  extract  pertinent  infor- 
mation such  as  line  amplitude,  linewidth,  line  position,  and  rms  noise. 


A Dewar  transfer  tube  was  inserted  into  the  NMR  crossed-coil  probe 
and  the  sample  was  located  in  the  Dewar.  Vapor  from  liquid  helium  (or 
liquid  nitrogen)  flowed  through  the  Dewar  to  cool  the  sample.  The 
liquid  helium  (or  nitrogen)  was  delivered  from  a storage  Dewar  by  an 
Air  Products  transfer  tube.  An  electronically  controlled  heater  raised 
the  temperature  of  the  vapor  to  the  desired  value  and  stabilized  it  by 
means  of  a temperature  sensor.  The  temperature  of  the  sample  was  moni- 
tored separately.  The  sensors  were  germanium  or  platinum  resistors  and 
gold-doped  iron  vs . chromel  or  copper  vs . constantan  thermocouples  de- 
pending on  the  particular  range  of  temperatures . 


Since  the  glass  Dewar  and  the  sample  were  both  located  in  the  NMR 
probe,  the  29Si  in  the  glass  also  contributed  to  the  resonance  line. 
The  magnitude  of  the  glass  contribution  to  magnetization  was  about  a 
factor  of  ten  smaller  than  the  sample  contribution  and  provided  a 
reference  sample  to  calibrate  the  spectrometer  gain. 


The  samples  were  procured  from  Heliotek,  Division  of  Textron, 
Sylmar,  California,  and  were  furnished  GFE.  They  were  prepared  from 
quartz  crucible  (Czochralski)  grown  silicon  containing  approximately 
1018  0/cm'3. 


Single  crystal  slabs  were  prepared  from  each  of  three  phosphorus 
doped  silicon  boules  with  0.1,  1,  and  10  fi-cm  room  temperature  resis- 
tivities, respectively.  The  slabs  have  the  dimensions  1 x 2 x 0.025  cm 
with  the  (lOO)  direction  approximately  normal  to  the  1 x 2 cm  surface. 
Representative  samples  from  each  of  the  "starting"  resistivity  cate- 
gories were  subsequently  lithium  diffused  (by  Heliotek)  to  concentra- 
tions from  1018  to  1019  Li/cm  % each  sub-category  containing  30  slabs. 


The  sample  resistivities  prior  to  lithium  diffusion  correspond  to 
phosphorus  concentrations  of  approximately  8.6  x 10*'  P/cm  (0.1  fi-cm) 
5.4  x 1018  P/cm’3  (l  fi-cm),  and  4.8  x 1014  P/cm'  (10  fi-cm).  These  con- 
centrations were  obtained  from  Hall-conductivity  techniques. 


The  samples  were  identified  by  two  numbers,  the  first  being  the 
starting  resistivity  and  the  second  the  logarithm  (base  10)  of  the 
lithium  concentration.  For  example,  the  sample  10-17  had  a starting 
resistivity  of  10  -cm  and  was  lithium  diffused  to  a concentration  of 
10 1 7 Li/cm  . 

Five  samples  were  selected  for  this  investigation;  i.e.,  0.1-15, 

10-17,  10-15,  1-15,  1-0. 

An  NMR  sample  consisted  of  fifteen  1 x 2 x 0.025  cm  single  crystals 
each  separated  by  a plastic  film  to  reduce  eddy  currents.  This  con- 
figuration was  necessary  to  avoid  skin  depth  problems  and  to  allow 
balancing  of  the  NMR  probe  for  higher  conductivity  samples . For 
measurements  taken  below  room  temperature,  a thermocouple  (chromel  vs. 
gold,  0.07  atomic  percent  iron)  was  securely  attached  to  the  sample. 

The  nuclear  spin-lattice  relaxation  time  for  each  sample  was 
measured  before  irradiation  in  the  temperature  range  from  8 to  300  K. 

The  samples  were  then  irradiated  at  Wright-Patterson  Air  Force  Base 
using  a High-Voltage  Engineering  Van  de  Graff  accelerator.  Pairs  of 
single  crystals  were  irradiated  (normal  to  the  1 x 2 cm  surface)  with 
1 MeV  electrons  and  were  maintained  at  room  temperature  during  the 
irradiation  using  compressed  air.  The  electron  beam  was  adjusted  to 
produce  uniform  F-center  coloration  of  a glass  plate  positioned  in  the 
sample  holder  and  the  fluence  level  was  determined  from  the  beam  cur- 
rent passing  through  the  sample.  Following  irradiation,  the  samples 
were  immediately  cooled  to  liquid  nitrogen  or  dry  ice  temperatures  to 
inhibit  lithium  annealing. 

The  saturation-recovery  method  was  used  for  Tj.  measurements . An 
adiabatic  fast  passage  to  the  center  of  the  resonance  (A  - 0)  aligns 
the  magnetization  with  the  rotating  (effective)  field.  Components  of 
the  magnetization  orthogonal  to  the  effective  field  decay  to  zero  in  a 
time  on  the  order  of  T ■ . Using  pulse  equipment  and  circuits  provided 
in  the  transmitter,  the  phase  of  the  rf  is  changed  by  rr/2  and  hence  the 
magnetization  is  orthogonal  to  the  effective  field.  Rf  phase  shifts  of 
it/2  rad  (shifting  time  t « T-,)  for  a time  greater  than  T^  will  lead  to 
saturation  of  the  ‘Si  system.  A series  of  fifteen  10  ms  rf  phase 
shifts  with  1.72  1 field  modulation  saturated  all  signals  observed  in 
the  equipment . A gate  signal  from  electronic  control  equipment  shifted 
the  static  field  to  the  center  of  the  line  and  turned  on  the  pulse 
generators.  After  saturation,  the  static  field  was  shifted  25  G below 
the  resonance  field  and  the  magnetization  was  allowed  to  grow  toward 

equilibrium.  line  transmitter  remained  on  during  this  recovery.  No  i 

difference  in  the  recovery  rate  was  observed  for  static  field  shifts  of 

25  or  50  G below  the  center  of  the  NMR  line.  After  a partial  recovery 

the  derivative  of  the  dispersion  signal  was  recorded  by  an  adiabatic 

fast  passage.  This  sequence  of  operations  including  saturation  was 

repeated  for  different  recovery  times . Figure  1 is  a typical  recording 

of  this  magnetization  recovery  information. 
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The  spin-lattice  relaxation  time  for  the  silicon  system  in  the  NMR 
probe  and  pyrex  sample  holder  dewar,  is  less  than  60  s for  all  tem- 
peratures and  frequencies  used  in  this  investigation.  The  magnitude  of 
this  background  signal  as  a function  of  temperature  was  measured  using 
experimental  conditions  similar  to  those  used  for  the  samples.  The 
curves  for  the  three  frequencies  wer  normalized  to  the  10.457  MHz 
background  signal  at  room  temperature  (45  arbitrary  units)  and  are 
exhibited  in  Fig.  2.  Both  the  8.464  and  10.457  MHz  data  were  obtained 
using  the  Varian  4-8  MHz  probe.  The  6.027  MHz  data  were  obtained  using 
the  Varian  2-4  MHz  probe.  The  majority  of  measurements  at  8.464  and 
10.457  MHz  were  made  using  rf  levels,  2HL , in  the  range  from  0.75  to 
1.25  G.  These  data  (Fig.  2)  were  used  to  normalize  the  sample  equi- 
librium magnetization  data  to  allow  its  temperature  dependence  to  be 
studied.  All  growth  times  used  in  obtaining  the  recovery  data  are  more 
than  four  times  TLg,  and  thus,  the  resulting  recovery  curve  appears 
approximately  as  an  exponential  offset  by  the  equilibrium  magnetization 

M_,  of  the  glass  background  ,9Si  system. 

& 

In  the  temperature  regions  where  the  recovery  exhibited  exponential 
behavior,  several  numerical  techniques  were  used  to  obtain  T]__1  from 
the  recovery  information.  Each  of  the  techniques  is  developed  from  the 
"ideal"  recovery  curve 

M(t)  = M0  - Ms  exp  (-t/Tj.)  = Mg  + Ms[l  - exp  (-t/Tj]  , (71) 

where  Ms  - Mo  - Mg.  The  quantities  Mg,  Ms,  and  Mq  are  the  equilibrium 
magnetizations  for  the  background,  sample,  and  combined  (background  and 
sample)  9Si  systems,  respectively. 

The  primary  data  reduction  technique  used  in  the  study  was  based 
on  a method  credited  to  Mangelsdorf . 01  By  using  equally  spaced  time 
intervals,  Eq.  (71)  can  be  rewritten  for  the  times  t = t0,  t0  + T,  ***, 


+ kT,  tr 


(k  + 1)T  , 


The  quantity  T is  the  interval  and  there 


is  no  restriction  on  the  initial  time  tG.  I 
tfc+i  t0  + (k  + l)T, 

M(tk)  = Mq  - Ms  exp  (-tk/Ti)  , 


Defining  tk  - t0  + kT  and 


M(tk  + T) 


Ms  exp  [-(tk+1)/Til  , or 


M(tk+1)  = exp  (-T/Ti)M(tk)  + Mq[1  - exp  (-T/Tl)]  . (72) 

A plot  of  M(tk+l ) vs.  M(tk)  will  give  Ti-1  from  the  slope; 

1/Tj.  -in  ( slope )/T  . (73) 

The  equilibrium  magnetization  Mo  is  obtained  from  the  intercept  of 
Eq.  (72)  with  the  equation  M(tk+L)  = M(tk).  Mg  is  obtained  by  extrapo- 
lating Eq.  (72)  to  M(0).  The  advantage  of  this  technique  is  that 
measurements  are  reduced  to  linear  graphs  and  standard  least-square 
fitting  techniques  can  be  applied. 
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Figure  2 . Background  ,3Si  signal  vs . temperature 


EXPERIMENTAL  RESULTS  AND  ANALYSES 


The  measured  values  of  TL  and  Ms  appear  in  Tables  1-11.  Values 
for  percent  error  were  obtained  from  the  rms  deviation  of  measured  M(t) 
from  the  best  fit  exponential  growth  curve  for  M(t).  Log-log  plots  of 
l/Ti  vs.  temperature  and  of  Ms  vs.  temperature  are  shown  in  Figs.  3-20. 

Figure  3 is  a graph  of  Ti-1  vs.  T for  samples  before  irradiation 
and  at  10.457  MHz.  Results  exhibited  in  Fig.  3 demonstrate  the  sensi- 
tivity of  the  nuclear  relaxation  rate  to  donor  impurity  concentrations . 
The  phosphorus  impurity  is  particularly  effective  in  relaxing  the  9Si 
system  at  low  temperatures.  This  is  indicated  by  an  increase  in  the 
slope  of  the  relaxation  rate  curves  when  the  phosphorus  concentration 
exceeds  the  lithium  concentration,  suggesting  that  samples  containing 
a substantial  concentration  of  lithium  tend  to  have  a greater  variety 
of  trap  energies  than  lithium- free  or  lithium-lean  samples. 
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The  results  for  the  sample  10-17  exhibit  a maximum  spin-lattice 
relaxation  rate  near  17  K.  The  absence  of  T^1  information  below  20  K 
for  all  other  samples  needs  consideration.  A loss  in  the  measured 
equilibrium  magnetization,  Ms,  accompanies  the  increasing  relaxation 
rates  with  decreasing  temperature.  This  loss  in  Ms  is  exhibited  in 
Fig.  4 for  the  corresponding  relaxation  rate  information  in  Fig.  3« 

The  background  ’’Si  signal  is  also  included  in  Fig.  4.  Several  obser- 
vations can  be  made  from  a comparison  of  these  figures; 

(1)  the  maximum  Ms  for  each  sample  occurs  near  the  maximum 
negative  slope  of  l/Tj_  vs.  T for  the  corresponding 
sample, 

(2)  the  samples  containing  higher  lithium  concentrations 
show  Ms  (max)  at  lower  temperatures , 

(3)  the  M,,  measurements  for  the  sample  10-17  indicate  a con- 
stant" or  minimum  equilibrium  magnetization  below  20  K, 
and 

(4)  the  background  °9Si  signal  is  less  temperature  dependent 
than  the  sample  signal  below  50  K. 

Attempts  to  obtain  Ms  and  Ti-1  for  samples  0.1-15  and  1-15  at  22  K 
and  20  K,  respectively,  were  unsuccessful.  The  precision  of  the  recovery 
information  used  in  determining  Ms  and  Ti-1  is  influenced  by  the  tem- 
perature dependences  of  these  quantities.  The  slope  of  the  Ti “ 1 vs.  T 
and  of  Ms  vs . T in  this  region  was  too  large  to  permit  accurate  measure- 
ment. At  lower  temperatures  no  NMR  line  could  be  displayed  except  for 
the  10-17  sample . 


Neither  the  increasing  relaxation  rates  nor  the  decreasing  equilib- 
rium magnetizations  with  decreasing  temperature  from  50  to  20  K have 
r r ■«=-.  ion-l  v recei  “d  attention  in  the  literat\ire  for  semiconducting 


i 

* 

K 


4 * 


?1Si  spin-lattice  relaxation  rate  vs.  temperature 
all  samples  before  irradiation,  10.457  MHz 


Figure  3 


Background 
29Si  Signal 


Figure  4.  Equilibrium  magnetization  vs.  temperature 
all  samples  before  irradiation,  10.457  MH 


vs.  T:  sample  0.1-15  before  irradiation 
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:,3i  spin-lattice  relaxation  rate  vs.  temperature 
1-0  before  irradiation 


Figure  7 


Vi  spin-lattice  relaxation  rate  vs.  temperature 
10-15  before  irradiation 


Figure  10 


spin-lattice  relaxation  rate  vs.  temperature 
before  irradiation 


Figure  11.  Equilibrium  magnetization  vs.  temperature 
10-17  before  irradiation 


vs.  T:  sample  10-1?  before  irradiation 
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nSi  spin-lattice  relaxation  rate  vs.  temperature 
0.1-15  after  irradiation,  before  annealing 


Figure  13 


:)Si  spin-lattice  relaxation  rate  vs.  temperature 
1-15  liter  irradiation,  before  annealing 


Fi  glare  15 


*Si  spin-lattice  relaxation  rate  vs.  temperature 
10-15  after  irradiation,  before  annealing 


^Si  spin-lattice  relaxation  rate  vs.  temperature 
10-17  after  irradiation,  before  annealing 


Before  Irradiation 
After  Irrcdiaticn 
After  Annealing 


TEMPERATURE (K) 


Integrate.!  electronic 
sample  10-17 


Figure  19 


" “Si  SDin-lattice  relaxation  rate  vs.  temperatur 
10-17  after  irradiation  and  13  hours  annealing 


silicon.  These  effects  have  been  observed  in  all  samples  considered  in 
the  present  investigation  and  are  sensitive  to  frequency  (or  field)  and 
paramagnetic  concentration.  Similar  behavior  has  been  observed  from 
the  °Si  background  signal  in  the  present  experiment . 

The  sample  0.1-15  is  doped  with  approximately  8.6  x lO1^  P/cm  ' and 
10*5  Li/cm  . For  this  sample  Hall-conductivity  measurements  and  EPR 
spectra  indicate  that  the  majority  of  donor  centers  are  ionized  above 
70  K and  the  phosphorus  system  provides  the  dominant  low  temperature 
paramagnetic  effects.  (The  conduction  electron  concentration  is  reduced 
by  an  order  of  magnitude  from  room  temperature  to  70  K.)  The  relaxation 
rate,  T^'1  vs.  temperature  results  for  the  three  operating  frequencies 
(and  fields)  used  in  this  research  are  exhibited  in  Fig.  5-  The  field 
dependence  and  the  increasing  relaxation  rates  with  increasing  tempera- 
ture between  50  and  300  K indicate  that  both  localized  and  nonlocalized 
electron  effects  are  contributing  to  the  relaxation  rate. 

Uncertainty  estimates  for  the  quantity  larger  than  5% 

are  indicated  by  error  bars  on  Fig.  5-  These  estimates  represent  the 
standard  deviation  in  the  T1"1  results  obtained  from  numerical  data 
reduction . 

The  nondegenerate  conduction  electron  contribution  [Eq.  (45)]  to 
T}.-1  is  field  independent  and  can  be  subtracted  from  the  experimental 
results.  Since  the  quantity  Ti”1  is  field  independent  at  room  tempera- 
ture, an  empirical  expression  'based  on  Eq.  (45)]  can  be  written  to 
describe  the  nondegenerate  conduction  electron  contribution  to  the 
relaxation  rate.  This  expression  is 

1/Tjc  = 1.35  x 10_?1NoT1  , (74) 

where  N0  is  the  mobile  carrier  concentration.  Equation  (74)  is  sketched 
in  Fig.  5 in  the  temperature  region  where  N0  has  been  determined  from 
Hall-conductivity  measurements  on  representative  single  crystals  from 
this  subcategory.  At  room  temperature  the  results  for  N agree  (to 
within  109'; ) with  the  concentrations  stated  by  the  sample's  source 
(Heliotek) . 

Several  authors  ”13  have  obtained  the  numerical  constant  in 
Eq.  (74)  from  room  temperature  data  and  concentrations  greater  than 
10!  mobile  carriers/cmJ ; their  results  give  the  numerical  constant 
1.09  x 10-  1 * 1 0$  Lampel  and  Solomon15  have  used  the  Overhauser 

effect  to  separate  the  conduction  and  localized  electron  contributions 
to  the  quantity  Ti"  . For  concentrations  between  10  and  101”  mobile 
carriers/cm  , they  found  the  conduction  electron  contribution  to  be  pro- 
portional to  N , at  room  temperature.  The  results  exhibited  in  Ref.  l4- 
16  indicate  that  localized  paramagnetic  contributions, 

1/Tj = l/Ti  - l/Tj.)c  to  the  quantity  T].-1  become  observable  for  mobile 
carrier  concentrations  smaller  than  10'  cm”'. 


i At  60,  90 5 135  and  195  K the  mobile  carrier  concentration  has  been 

obtained  from  Hall-conductivity  data.  These  data  and  Eq.  (74)  allow  an 
estimate  of  the  localized  electron  contribution  to  the  relaxation  rate 
to  be  made.  In  Fig.  6,  the  localized  electron  contribution  (Til-1)  to 
the  relaxation  rate  is  plotted  as  a function  of  temperature  between 
60  and  200  K.  The  broad  maxima  in  these  curves  indicate  that  Tg  is  not 
a strong  function  of  temperature.  Further,  the  maxima  in  Fig.  6 indi- 
cate that  tc  is  of  the  same  order  of  magnitude  as  l/up  near  135  K 

A field-dependent  contribution  to  the  quantity  Ti-1  is  described 
by  Eqs.  (49),  (53)  5 (59)  > and  (64).  The  rapid  diffusion  case  is 
expected,  since  Ti  » To  [see  Eqs.  (64)  and  (66)].  These  equations  are 
summarized  in  Eq.  (75)  and  describe  relaxation  by  localized  paramagnetic 
centers ; 

j 1/T,.l  « (8Jr/l5)(un/<uet>)3/4«S'yi7l[S(S  + l)/a  ]Np[Tc/(l  + u£t|)] 


« 4.03  x 10" 1 1 (t/Hg )3/4Np[ tc/ (l  + wgr |)]  , (75) 

where  Np  is  the  localized  paramagnetic  concentration. 

At  8.464  MHz  (H  , - 10.0  kG)  and  135  K,  Ti"1  = 5-65  x 10~4  s'1. 

Eq.  (75)  then  gives 
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Np(Tc/(1  + wLTc) ) ^ 3-54  x 10"  s-cm"'  . (76) 

The  quantity  ic/(l  + w£t|)  has  a maximum  (ic/2)  when  tc  - l/wp.  Thus, 
the  minimum  localized  paramagnetic  concentration  predicted  by  Eq.  (76) 
for  this  temperature  is  3-8  x 101  cm*  '.  The  contribution  of  these 
localized  centers,  Til"1,  3s  of  the  same  order  of  magnitude  as  the  non- 
degenerate conduction  electron  contribution  to  Tj.'1.  The  number  of 
electrons  localized  to  the  donor  impurity  system  at  135  K is  approxi- 
mately 4.9  x 10:’’  cm'  and  represents  roughly  half  of  the  donor  im- 
purities ionized. 

To  apply  the  rapid  diffusion  model,  the  quantity  5 - p' /2b  « 1. 

The  scattering  length  p - (C/D)l/,  ~ 6. 91  x 10'  cm,  where  the  nuclear 
spin-diffusion  constant  (D  - 2 x 10' • cm  s"1)  was  obtained  from  the 
work  of  Jerome  and  Winter.  The  correlation  time  used  to  calculate  C 
[Eq.  (53)  ]»  was  taken  as  -tc  1.88  x 10"  s.  The  diffusion  barrier 
radius,  b 2.73  x 10"  cm  was  obtained  from  Eqs.  (59)  and  (49).  These 
numerical  results  yield  0.03  for  the  quantity  5.  The  field  dependence 
in  Fig.  6 is  in  qualitative  agreement  with  Eq.  (75)-  Between  90  and 
60  K and  at  10.457  MHz,  the  increasing  relaxation  rates  with  decreasing 
temperature  indicates  that  a paramagnetic  system  is  becoming  "active" 
in  nuclear  re] axation  below  70  K. 

The  linewidth  barrier  radius,  l , for  the  localized  paramagnetic 
system  can  be  calculated  from  Eqs.  (49)  and  (6l).  Since  tc  x l/io^,, 
t 4.5  x 10  (Hq/T)1*  ox  20  ' at  13  and  with  a 12.4  kG  • a tic 


to 
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field.  This  result  represents  approximately  0.03%  of  the  3Si  nuclei 
being  removed  from  the  NMR  line  at  135  K.  The  observed  equilibrium 
magnetization  for  the  ~*Si  system  should  thus  increase  according  to  the 
Curie  law  (T*1).  The  curves  in  Fig.  4 indicate  an  approximate  Curie  law 
behavior  between  50  and  300  K. 

The  temperature  dependence  of  tc  (indicated  by  the  broad  maxima 
in  Fig.  6)  is  not  as  strong  as  would  be  expected  for  a Raman  type 
electron  spin-lattice  relaxation  mechanism. 9P  Further,  by  comparing 
Figs.  3 and  6,  the  quantity  Tip,”1  for  the  sample  0.1-15  is  nearly  three 
times  larger  than  the  corresponding  relaxation  rates  (Tx-1)  for  the 
samples  containing  approximately  5 x 1015  donor  impurities/cmJ  These 
results  imply  that  the  electron  concentration  localized  to  the  donor 
impurity  system  is  contributing  to  the  nuclear  relaxation  between  50 
and  300  K . 

Below  50  K,  the  majority  of  conduction  electrons  are  localized  to 
the  donor  impurity  system.  The  EER  spectra  of  this  sample  reveal  that 
the  phosphorus  donor  impurity  is  the  dominant  paramagnetic  system  in 
this  temperature  range.  The  increasing  relaxation  rates  and  the  loss 
in  the  equilibrium  magnetization  with  decreasing  temperature  are  the 
gross  NMR  effects  introduced  by  the  phosphorus  donor  system  (see 
Figs.  5 and  4). 

An  auxiliary  experiment  was  performed  on  sample  0.1-15  at  6.027  MHz 
involving  changes  in  temperature  between  32  K and  6 K in  times  short 
compared  with  T;  for  the  sample.  First  Tx  and  Ms  were  measured  at  32  K 
and  it  was  noted  that  Ms  at  32  K was  in  the  Curie-law  region  and  Tx  was 
~ 1000  s . Then  the  sample  was  saturated  at  32  K and  allowed  to  grow 
(at  32  K)  for  a time  (4000  s)  long  compared  with  Tx  to  achieve  M = Ms 
(32  K).  The  temperature  was  then  lowered  quickly  to  6 K and  allowed 
to  attain  thermal  equilibrium  (615  s)  at  that  temperature  at  which  time 
the  line  was  recorded.  The  resulting  line  amplitude  was  smaller  by  a 
factor  of  two  than  Ms  (32  K).  The  6 K recording  appeared  to  be  narrowed 
slightly  (10$)  and  shifted  40  mG  to  higher  fields  which  is  not  signifi- 
cant, however,  it  definitely  did  not  appear  to  be  broadened. 

As  a second  step  the  sample  was  saturated  at  6 K and  allowed  to 
grow  for  4615  s and  was  recorded.  A very  small  amplitude,  - 5%  of  Ms 
(32  K)  was  obtained.  This  can  be  used,  along  with  the  previous  data, 
to  estimate  Tx  at  6 K to  be  ~ 80  hours  or  about  288  times  longer  than 
Tx  (32  K),  although,  the  estimate  is  questionable  since  it  involves  an 
estimate  of  Ms  (6  K)  by  use  of  the  Curie-law  relation.  It  is  evident, 
however,  that  TL  (6  K)  is  very  long  compared  with  Tx  (32  K). 

In  the  third  step,  the  sample  was  saturated  at  6 K,  allowed  to 
grow  at  bK  for  4000  s,  raised  quickly  in  temperature  to  32  K and 
recorded  after  615  s (4615  s total)  and  was  found  to  have  an  amplitude 
consistent  with  almost  615  s growth  time  at  32  K,  based  upon  Ms  (32  K) 
and  Tx  (32  K)  previously  measured. 
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On  the  basis  of  the  results  of  this  experiment  an  estimate  of  the 
linewidth  barrier  radius  can  be  made.  Since  approximately  l/2  of  the 
^Si  spins  were  removed  from  the  resonance  line  after  cooling  from  32  K, 
the  linewidth  barrier  radius  can  be  obtained  from  the  paramagnetic  con- 
centration. For  sample  0.1-15  it  is  approximately  8.7  x 10 1 impurities/ 

cm  . Thus,  i R/'^2  (l40  X/V I)  111  X.  This  result  can  be  com- 

pared with  0 obtained  using  Eq.  (6l)  and  the  static  electron  magnetic 
moment ; 


l * a[ue/un]1/3  (77) 

where  the  quantity  a is  the  'Si  separation; 

a 7.43  X , 

ae  9*27  x 10-i?1  erg/G  , and 

un  = 2.80  x lO-'"’4  erg/G  . 

Equation  (77)  then  gives  0 111  X>  which  is  in  excellent  agreement  with 

the  results  obtained  from  the  impurity  concentration.  These  results, 
when  combined  with  the  Tj._!  and  the  Ms  information  obtained  between 
20  and  30  K,  indicate  the  role  of  the  localized  electron  relaxation 
rates  T^e-  and  T^g-1  in  the  nuclear  relaxation.  The  electron  spin- 
spin  relaxation  rate  for  the  localized  electron  system  will  provide  an 
upper  bound  on  the  electronic  correlation  time.  An  estimate  of  T-e  is 
obtained  from  the  electron  magnetic  dipole-dipole  interaction  and  the 
impurity  concentration; 

l/T2e  - 7ehLe  2p§/fcR  or  Toe  0.14  ms  . (78) 


The  electronic  correlation  time  is  of  the  came  order  of  magnitude  as 
the  nuclear  spin-spin  relaxation  time; 

l/Tp  - 7 h (5.32  x 103  md/Gs) (0.5  G)  or  T?  =*  O.38  ms  . 

(79) 

For  a ve  5 ety  of  impurities  including  phosphorus,  Castner  was 
able  to  fit  his  electron  spin-lattice  relaxation  data  to  an  empirical 
expression  of  the  form 


1/T 


1 e 


AHqT  + BH  t'  + CT1  + DT13  + E(H)  exp  (-A/AT) 


(80) 


where  the  term  AH' T represents  a direct  phonon  (emission  or  absorption'! 
process  and  the  remaining  t'-rms  represent  Raman  type  processes.  These 
mechanisms  a^cr  discussed  in  some  detail  in  the  book  by  Standley  and 
Vaughan  and  their  references.-1  Castner  has  measured  Tle-1  in  Czochralski 
grown  silicon  with  1.5  x 10  14  l/cm ' (1.5  fi-cm)  for  temperatures  in  the 
range  from  1.5  to  30  K.  (This  sample  also  contained  2.5  x 1013  As /cm.) 
Above  6 K his  results  fit  the  empirical  expression 


LV-T1 
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1/Tj.e  = E exp  (-A/kT) 


(81) 


where  E 2.2  x 10^  s'1  and  A 11.6  MeV.  The  quantity  A is  the 
activation  energy  for  the  IS  singlet  ground  state  for  the  phosphorus 
donor  impurity.  The  activation  energy  was  found  to  be  the  same  in  a 
Merck  float  zoned  silicon  crystal  with  9 x 101*5  P/cm'  (0.7  ! -cm). 

Castner’s  measurements  of  the  electron  spin-lattice  relaxation 
rate  indicated  transition  times  of  the  order  of  millisecond:  near  20  K. 
'Thus,  the  loss  in  equilibrium  magnetization  exhibited  in  Fig.  ■ may  be 
associated  with  the  growth  of  the  linewidth  barrier  radius.  Further, 
the  long  nuclear  relaxation  time  observed  at  6 K may  be  1 he  result  of 
correlation  times  for  the  local  field  cf  the  ordei  f the  r .el«n  r spin- 
spin  relaxation  time. 

The  increasing  relaxation  rate  with  decreasing  te:.v.  erature 
(20  K < T < 50  K)  for  sample  0.1-15  can  be  describe  1 fc  the  em  i.rical 
expression 

I/Tj.  - 5.87  x lO-15  exp  (82.6/t)  . (82) 

Equation  (82)  is  included  in  Fig.  5 and  implies  that  the  average  corre- 
lation time  for  the  local  field,  t , is  proportional  to  the  electron 
spin-lattice  relaxation  time,  Tie-*  Since  the  observed  relaxation  rate 
(Tx-')  increases  with  tc  and  is  field  independent,  then  tc  < l/up.  The 
reduction  of  the  magnetization  with  decreasing  temperature  and  the  in- 
creasing relaxation  rate  with  decreasing  temperature  imply  that  a dis- 
tribution of  correlation  times  about  the  average  tc  is  responsible  for 
the  NMR  behavior  between  20  and  30  K;  that  is,  tc  is  a function  of  both 
position  in  the  sample  and  temperature  of  the  sample.  Nuclear  spins 
with  longer  tc  are  lost  from  the  NMR  line,  while  spins  which  experience 
a local  field  with  shorter  correlation  times  contribute  to  the  observed 
magnetization . 

The  samples  1-0,  1-15,  and  10-15  each  contain  approximately  5 x IG15 
donor  impurities/cm".  As  with  the  sample  0.1-15,  the  EPR  spectra  of 
these  samples  demonstrate  the  presence  of  localized  paramagnetic  sites 
below  50  K.  The  NMR  relaxation  rate  results  are  exhibited  in  Figs.  7» 

8,  and  9. 

Above  50  K,  the  relaxation  rates  for  these  samples  are  nearly  an 
order  of  magnitude  smaller  than  those  for  the  sample  0.1-15 . As  a 
result  of  the  lower  donor  impurity  concentration,  the  frequency-dependent 
contribution  (above  50  K)  is  resolved  at  room  temperature.  At  60  K,  the 
field-dependent  relaxation  rate  is  nearly  an  order  of  magnitude  larger 
for  the  sample  0.1-15  than  for  the  samples  considered  in  this  section. 
This  result  indicates  that  the  correlation  time  decreases  and/or  the 
localized  paramagnetic  concentration  increases  with  increasing  donor 
impurity  concentration.  A combination  of  these  effects  is  expected. 
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The  room  temperature  mobile  carrier  concentration  for  the  sample 
1-0  has  been  obtained  from  Hall-conductivity  data  and  is  near  4.8  x 10i8 
electrons/cm3.  From  this  information  and  Eq.  (74),  the  nondegenerate 
conduction  electron  contribution  to  Tj.~‘  is  approximately 
l/Ti)  -•«  1.1  x 10' 1 s'1.  By  subtracting  this  contribution  from  the 
measured  room  temperature  relaxation  rates  for  the  sample  1-0,  field- 
dependent  contributions  remain  which  are  of  uhe  same  order  of  magnitude 
as  the  experimental  error  at  room  temperature  for  the  sample  0.1-15. 

Near  room  temperature  for  the  sample  0.1-15  the  number  of  electrons, 
localized  to  the  donor  impurity  system,  is  roughly  5 x 101"  e'/cm '. 

Heat  treatment,  accompanying  lithium  diffusion,  produced  a two-fold 
increase  in  the  relaxation  rates  for  the  sample  1-15  over  those  for  the 
sample  1-0. 

At  room  temperature  for  concentrations  smaller  than  10 37  donor 
impurities/cm  % several  authors  have  demonstrated  larger  relaxation 
rates  than  could  be  predicted  by  the  nondegenerate  conduction  electron 
theory  [Eq.  (45)]. 14-16  The  elemental  donor  impurity  and  sample  pre- 
paration (including  oxygen  content)  were  found  to  have  little  effect  on 
the  low  concentration  relaxation  rates.  Rahilly,16  observed  a field 
dependence  at  room  temperature  in  Czochralski  grown  samples  with  fewer 
than  101  donor  impurities/cm ' . (After  subtracting  the  nondegenerate 
conduction  electron  contribution,  Tx  was  found  to  be  proportional  to 
a)' ) In  a sample  containing  10JS  Li/cm  (near  the  semiconductor  to 
metal  transition),  Rahilly  found  the  quantity  TL-1  independent  of  field 
at  room  temperature  and  then  to  become  field  dependent  at  lowered  tem- 
peratures. The  field-dependent  relaxation  rates  between  20  and  300  K 
for  samples  with  donor  impurity  concentrations  less  than  1018  cm-3  have 
not  previously  been  reported  in  the  literature  for  semiconducting 
silicon. 

The  field  dependent  results  for  these  samples  indicate  that  the 
time  dependence  of  the  local  field  may  be  influenced  by  the  interaction 
of  conduction  electrons  with  the  donor  impurity  sites . These  inter- 
actions produce  power  spectra  for  the  local  field  of  the  form 

j(u)  ~ XbF t|  { [ sin  (unb/2 ) ]/(coxb/2 ) (83) 
or 

j(w)  - b T j/ (l  + CO  Tj)  , (84) 

where  xb  is  the  residence  time  for  an  electron  in  a trap  site,  xj  is 
the  carrier  lifetime,  and  A is  the  trapping  frequency  at  a site.  If  xb 
and/or  x-i  are  short  compared  to  l/co  , the  power  spectrum  will  be  fairly 
uniform  (white).  This  is  the  high  temperature  and  concentration  limit. 
Then,  the  quantity  Tj.-1  will  be  field  independent  and  determined  by  b 
or  1 ;,  and  b^.  When  Tb,  Tj  > l/wb,  a field  dependence  will  be  observed 
which  contains  information  on  the  power  spectrum  of  the  local  field. 

The  quantity  A is  controlled  by  the  kinetics  of  the  conduction  electon 
system.  The  quantities  tb  and  t.j  depend  on  temperature,  trap  depth. 
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nd  lossitly  field.  Shallow  donor  states  have  wave  functions  extending 
over  many  1 ttice  constants,  and  thus,  the  carrier  lifetime  is  very 
short  for  these  states.  Typical  carrier  lifetimes  are  on  the  order  of 
microseconds,1  and  hence,  modulations  in  the  local  field  will  have 
frequency  components  near  the  IIMR  angular  frequency.  Finally,  the 
quantity  b represents  the  rms  amplitude  of  the  local  field  modulations 
due  to  carrier  trapping  at  donor  sites.  The  spectra  described  by 
Eqs.  (83)  and  (84)  are  similar  to  that  for  localized  paramagnetic  elec- 
trons 'Eq.  (75)1* 

As  with  the  sample  0.1-15,  localized  paramagnetic  effects  become 
important  below  50  K.  The  sample  1-0  was  not  doped  with  lithium,  but 
contains  5*4  x 1015  p/cm*.  This  sample  exhibits  the  largest  temperature 
dependences  in  the  relaxation  rate  below  50  K.  The  dashed  curve  in 
Fig.  7 follows  the  empirical  expression 

1/Tl  1.02  x 10-"  exp  (134/t)  . (85) 

The  constant  in  the  exponential  (134  K or  12.7  MeV)  is  approximately 
the  same  as  Castner’s  result9*  for  the  excitation  energy  for  the  ground 
state  of  the  phosphorus  donor  impurity  (122.5  K or  11.6  MeV).  His  re- 
sults were  obtained  from  electron  spin-lattice  relaxation  rate  measure- 
ments . 

Below  30  K,  two  measurements  of  TL_i  for  the  sample  1-0  at 
8.464  MHz  indicate  a deviation  from  the  empirical  expression  (85).  A 
similar  datum  is  reported  for  the  sample  0.1-15  at  25  K.  These  results 
describe  the  long-term  recovery  of  the  magnetization  in  a region  where 
the  recovery  is  not  exponential.  The  temperature  at  which  the  non- 
exponential recovery  appears  (=  28  K for  samples  0.1-15  and  1-0 ) is 
approximately  the  temperature  where  Castner’s  results  indicate 
TLe  =»  l/wL*  The  short-term  recovery  could  not  be  separated  from  the 
experimental  data  for  these  samples  in  the  28  K region. 

The  lithium  concentration  is  comparable  to  or  larger  than  the 
phosphorus  concentration  in  the  samples  1-15  and  10-15 . The  field 
dependence  from  the  conduction  electron  relaxation  region  continues 
below  50  K.  Also,  the  slope  of  the  relaxation  rate  curves  is  smaller 
for  these  samples  than  for  the  more  heavily  phosphorus  doped  samples . 
These  results  are  in  qualitative  agreement  with  the  25$  lower  ioniza- 
tion energy  of  the  lithium  donor  center. 

Further  evidence  for  the  involvement  of  the  lithium  donor  center 
in  the  nuclear  relaxation  is  provided  by  the  equilibrium  magnetization 
measurements  (Fig.  8)  Samples  with  higher  lithium-to-phosphorus  ratios 
approximate  Curie-law  behavior  to  lower  temperatures.  This  observation 
is  consistent  with  the  shorter  correlation  times  associated  with  the 
lithium  donor  system. 

The  most  interesting  results  for  the  sample  10-17  are  the  broad 
maxima  in  the  relaxation  rate  curves  below  20  K and  the  accompanying 
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low-temperature  behavior  of  the  equilibrium  magnetization.  These 
effects  are  exhibited  in  Figs.  10  and  11. 

The  nondegenerate  conduction  electron  (field  independent)  contri- 
bution was  obtained  from  Eq.  (74)  using  the  conduction  electron  concen- 
trations obtained  from  Hall  data.  This  contribution  is  sketched  in 
Fig.  10  and  has  been  removed  from  the  experimental  results  in  Fig.  12. 
The  maxima  for  these  curves  are  listed  below  along  with  an  estimate  of 
the  correlation  time,  t*  = l/ooj , at  these  temperatures. 


T 

T*max 

(s) 

M 

(s  x) 

(see  text) 

95 

1.55  x 10-:i 

2.64  x 10 “8 

105 

1.09  x 10-3 

1.88  x 10-8 

115 

6.6  x 10- ‘ 

1.52  x 10-8 

The  correlation  times,  t*,  represent  maxima  for  the  quantity 
xc/(l  + Tq ) and  not  maxima  for  the  quantity  Til-  as  given  by  Eq.  (75 ) - 
However,  the  broad  maxima  in  Fig.  12  indicate  that  tc  is  not  a strong 
function  of  temperature  and  the  quantity  may  be  used  for  order  of 
magnitude  calculations.  The  maxima  of  Fig.  12  occur  at  lower  tempera- 
tures than  those  for  the  sample  0.1-15.  This  result  is  consistent  with 
a decreasing  correlation  time  with  increasing  impurity  concentration. 

At  105  K the  number  of  electrons  localized  to  the  donor  impurity 
system  is  approximately  8.4  x 10'  cm-'.  This  result  can  be  compared 

with  the  paramagnetic  concentration  estimated  by  Eq.  (75); 

Np  = 8.8  x 10 15  cm-'. 

Hall -conductivity  measurements  and  EFR  spectra  on  sample  10-17 
indicate  that  the  majority  of  conduction  electrons  have  become  localized 
electrons  below  100  K.  The  effects  of  this  localized  system  on  the  KMR 
results  are  seen  below  40  K in  Figs.  10  and  11.  The  maximum  relaxation 
rates  below  20  K and  the  T-1  •3‘  (~  Curie  law)  dependence  in  the  equi- 
librium magnetization  indicate  that  the  average  correlation  time,  :c, 
is  of  the  order  of  0.02  |is  at  15  K.  Information,  obtained  from  the 
maxima  for  the  curves  in  Fig.  10  (below  20  K),  is  exhibited  below. 

(Data  below  20  K) 


T 

T1-1 

^c 

(K) 

(s-1)  x 10-8 

( s ) x 10- 

10.5 

2.64 

2.49 

13.5 

1.88 

1-93 

17.0 

1.52 

1 .36 
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These  results  can  be  fitted  to  the  empirical  expression 

1/Tl  C'[  ic/(1  + c^t®)]  , (86) 

where  C'  1.88  x 10s  s"1  ± 10%.  By  comparing  Eqs.  (75)  and  (86),  an 
estimate  of  the  localized  electron  concentration  can  be  obtained; 

Np  = 6.6  x 10* 7 cm-  . This  result  compares  with  1.1  x 101  cm'  for  the 
donor  impurity  concentration.  The  temperature  dependence  for  ic  above 
is  smaller  than  that  for  a direct  or  Raman  electron  spin-lattice  re- 
laxation process  and  may  be  the  result  of  a broad  distribution  of  cor- 
relation times  in  the  sample.  A distribution  of  correlation  times  about 
the  average,  ic,  is  expected  in  sample  10-17,  because  of  the  sample's 
nonuniform  lithium  concentration1’  and  the  tendency  of  lithium  to  form 
clusters  in  silicon.  14 

A distribution  of  correlation  times  is  also  implied  by  the 
equilibrium  magnetization  measurements . The  quenching  of  the  observed 
equilibrium  magnetization,  Ms,  is  not  complete  for  sample  10-17  (see 
Fig.  11  for  the  constant  or  minimum  equilibrium  magnetization  below  20  K). 
At  approximately  23  K the  local  field  components  with  longer  correlation 
times  begin  to  remove  nuclear  spin  contributions  from  the  NMR  line.  The 
ratio  of  the  temperatures  at  which  this  quenching  becomes  noticeable  for 
the  samples  0.1-15  and  10-17  is  0.71.  This  ratio  is  in  good  agreement 
with  the  ratio  of  the  ground-state  ionization  energies  for  the  donor 
centers  [Li-0  and  P],  0.75;  Raman  type  processes  are  thus  implied  for 
temperatures  above  20  K. 

For  -rc  as  1 /cop,  the  linewidth  barrier  radius  Z - 4.5  x 10“e  (Hq/t)" 
and  only  4$  of  the  3Si  spins  would  be  removed  from  the  NMR  line; 

(H0  - 10.0  kG,  T = 15  K,  and  Np  6.6  x 10 ' cm"').  However,  if  the 
transition  rate  for  some  of  the  electronic  spins  is  near  T2-1,  the 
total  electron  moment,  uez,  will  be  "seen"  by  the  ’Si  nuclei  and  the 
linewidth  barrier  radius  will  be  111  A.  The  results  in  Fig.  4 indicate 
that  approximately  50%  of  the  3Si  spins  are  removed  from  the  NMR  line 
by  15  K. 

The  constant  or  increasing  equilibrium  magnetization  observed 
below  20  K is  a consequence  of  a loss  in  the  integrated  electronic  para- 
magnetism observed  in  sample  10-17-  The  loss  in  paramagnetism  below 
20  K (see  Fig.  17)  is  attributed  to  the  formation  of  donor  pairs,  which 
become  diamagnetic  at  lower  temperatures.  The  excitation  of  these 
donor  pairs  (or  clusters)  to  the  paramagnetic  states  provides  spectral 
components  to  the  local  field. 

Many  of  the  changes  introduced  by  electron  irradiation  are  the 
result  of  carrier  removal  from  the  conduction  band  and  the  formation  of 
new  paramagnetic  sites.  These  changes  are  reviewed  in  terms  of  the 
preceding  discussions  and  existing  models  for  the  irradiation  damage. 

The  sample  0.1-15  was  irradiated  at  room  temperature  with  1 MeV 
electrons  to  a fluence  of  1.5  x 101’  e"/cm  . At  300  K the  * ’Si  relax- 
ation rate  was  unchanged  within  the  precision  of  the  measurements  (±  5%) • 
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The  characteristic  properties  of  the  phosphorus  donor  system  are  still 
present  in  the  relaxation  behavior  below  room  temperature  (see  Fig.  13). 

At  100  K and  10.457  MHz,  the  quantity  Tj.-1  is  approximately  60% 
smaller  than  T1-1  observed  before  irradiation.  The  difference  is  negli- 
gible at  6.027  MHz  at  100  K.  The  smaller  relaxation  rates  near  this 
temperature  indicate  the  trapping  of  conduction  electrons  at  a higher 
temperature. 

Below  70  K,  localized  paramagnetic  electrons  affect  the  nuclear 
relaxation  and  the  observed  equilibrium  magnetization.  Within  the 
precision  of  the  measurements  (as  large  as  + 35%)  the  results  at  8.464 
and  10.457  MHz  fit  the  empirical  expression  Eq.  (102)  for  the  relaxation 
rates  before  irradiation.  However,  the  6.027  MHz  datum  at  24  K is  50% 
larger  than  the  corresponding  T].-1  before  irradiation.  The  measurements 
at  24  K following  irradiation  provided  an  excellent  fit  (±  3%)  to  the 
recovery  information.  An  empirical  expression  fitting  the  experimental 
results  between  24  and  50  K is  written 

l/Ti  = 6.24  x 10"5  exp  (90.9/T)  . (87) 

The  average  binding  energy  indicated  by  this  expression  is  107  larger 
than  before  irradiation.  These  results  are  consistent  with  the  equi- 
librium magnetization  measurements  exhibited  in  Tables  1 and  6.  Before 
irradiation,  the  maximum  Ms  at  8.464  MHz  occurred  at  32  K and  following 
irradiation  (and  before  annealing)  the  maximum  occurred  at  approximately 
40  K. 


Samples  1-0,  1-15,  and  10-15  and  sample  10-17  were  irradiated  with 
1 MeV  electrons  to  a fluence  level  of  10  e“/cm  . At  65  K and 
6.027  MHz,  the  quantity  Tv  for  sample  1-0  was  found  to  be  17,000  s 
following  irradiation;  this  compares  to  approximately  12,000  s for  the 
sample  before  irradiation. 

The  post-irradiation  results  before  annealing  for  sample  1-15  are 
exhibited  in  Fig.  22.  A two-fold  increase  in  the  relaxation  rates 
(below  70  K at  6.027  MHz)  is  revealed  by  comparing  Figs.  16  and  22. 

The  increased  relaxation  rates  at  these  temperatures  are  an  indication 
of  deep  paramagnetic  trap  sites. 

Watkins  and  associates9'’  > 9 ' have  observed  EPR  spectra  resulting 
from  irradiation-induced  damage  sites  in  high-oxygen-content  phosphorus- 
doped  silicon.  Their  results  indicated  that  the  principal  paramagnetic 
defect  introduced  by  1 MeV  electron  irradiation  is  the  oxygen-vacancy 
pair  [A-center).  The  EFR  spectrum  for  this  center  has  been  observed  in 
these  samples. 

In  the  sample  with  the  smallest  controlled  dopant  concentration 
(10-15),  the  electron  irradiation  changes  the  NMR  relaxation  behavior 
at  all  temperatures.  These  changes  include  (l)  an  increase  in  the  re- 
laxation rates  by  approximately  4o%  at  200  K,  (2)  a decrease  in  the 


frequency  dependence  above  100  K,  (3)  a large  frequency  dependence  in 
the  relaxation  rates  between  6.02?  and  8.464  MHz,  and  (4)  a small  fre- 
quency dependence  between  8.464  and  10.457  MHz.  The  frequency  depend- 
ence appearing  in  Fig.  15  implies  a minimum  or  constant  power  spectrum 
in  the  8-11  MHz  frequency  range.  A similar  field  dependence  is  seen 
in  Fig.  9 near  room  temperature. 

Sample  10-17  was  given  greater  attention  following  electron 
irradiation  because  of  its  interesting  pre-irradiation  behavior.  The 
NMR  relaxation  results  are  exhibited  in  Fig.  16.  Following  electron 
irradiation,  the  relaxation  rates  near  room  temperature  are  reduced  by 
15$  and  the  frequency  dependence  below  20  K becomes  smaller. 

The  frequency  dependence  above  40  K extends  to  room  temperature, 
indicating  an  increased  localized  paramagnetic  concentration.  Since 
sample  0.1-15  did  not  indicate  a change  in  the  conduction  electron  re- 
laxation contribution  near  room  temperature,  the  lithium  donor  system 
may  be  responsible  for  the  changes  observed  in  sample  10-17.  The 
reduction  in  relaxation  rates  above  40  K following  irradiation  is 
attributed  to  a reduction  in  the  effectiveness  of  lithium  donor  elec- 
trons as  relaxing  agents,  possibly  because  of  a change  in  the  power 
spectrum  toward  shorter  correlation  times . 

Figure  17  is  a graph  of  the  integrated  electronic  paramagnetism 
(arbitrary  units)  vs.  temperature  for  sample  10-17.  Nearly  an  order  of 
magnitude  increase  was  observed  in  the  electronic  paramagnetism  at  25  K 
following  irradiation.  By  comparing  Figs.  16  and  17,  the  localized 
paramagnetic  system  is  seen  to  become  important  in  relaxing  the  nuclear 
system  below  40  K.  The  large  increase  in  electronic  paramagnetism  is 
not  reflected  in  the  NMR  Tj.-1  measurements,  possibly  because  of  an 
unfavorable  distribution  in  correlation  times . 

Brucker  and  associates9’-  have  proposed  a model  which  is  in  quali- 
tative agreement  with  the  NMR  measurements  presented  here  and  the  EPR 
measurements.  The  model  includes  changes  in  the  concentration  of  three 
major  centers  in  lithium-doped  quartz-crucible-grown  silicon  during 
electron  irradiation.  These  are 

(1)  the  lithium  Ll-center  described  by  the  equation, 

Li+  + V=  LiV"  ; 

(2)  the  lithium  L2-center  described  by  the  equation, 

Li0+  + V=  e LiOV"  ; and 

(3)  the  oxygen-vacancy  center  (A-center)  described  by  the 
equation , 

0 + V~  ^ 0V"  [A-center]  + e~ 
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The  quantities  V"  and  V~  are  singly  and  doubly  ionized  vacancy  state-. 
These  equations  indicate  the  formation  of  paramagnetic  centers  and 
changes  in  trap  energies  of  lithium  as  a donor  impurity. 

Following  electron  irradiation,  samples  0.1-15  and  10-17  were 
•annealed.  The  NMR  results  still  possess  the  gross  properties  of  a donor 
impurity  system.  However,  changes  in  the  relaxation  behavior  were 
observed  after  annealing  and  these  are  attributed  to  changes  in  the 
paramagnetic  defect  system. 

The  sample  0.1-15  was  annealed  for  60  hours  at  room  temperature 
and  the  NMR  results  following  annealing  are  exhibited  in  Fig.  l8. 

Below  50  K the  relaxation  rates  at  8.464  and  10.457  MHz  can  be  fitted 
to  the  Eq.  (82)  to  the  precision  of  the  data.  Above  40  K measurements 
of  the  quantity  TL_1  indicate  a smaller  field  dependence  between  8.464 
and  10.457  MHz  than  between  6.027  and  8.464  MHz.  This  was  not  observed 
before  annealing. 

The  effects  of  annealing,  demonstrated  by  sample  10-17,  are 
particularly  interesting  in  terns  of  the  irradiation  damage  and  lithium 
annealing  model  proposed  by  Brucker.34  Figure  19  is  a plot  of  the  NMR 
spin-lattice  relaxation  rate  results  following  a 13-hour  anneal.  The 
field  dependence  continues  below  30  K and  indicates  that  correlation 
times  longer  than  l/wL  are  affecting  the  relaxation  at  all  temperatures 
considered.  The  equilibrium  magnetization  data  at  each  frequency  before 
irradiation  and  after  irradiation  (before  and  after  annealing)  are 
exhibited  in  Fig.  20  for  sample  10-17- 

The  relaxation  rates  above  40  K are  reduced  further  (by  approxi- 
mately 40$  at  80  K)  following  annealing.  The  results  indicate  a 
change  in  the  effectiveness  of  lithium-related  donor  electrons  in  con- 
tributing to  relaxation. 

Brucker  and  associates  proposed  that  the  oxygen-vacancy  and  the 
lithium  LI  and  L2  centers  are  present  in  quartz-crucible-grown  lithium- 
doped  silicon  immediately  following  room  temperature  electron  irradi- 
ation. During  annealing,  lithium  forms  further  complexes  with  defect 
centers.  This  mechanism  is  described  by  the  equation 

Li+  + Li+  + V“[A-center]  + e-  — ne"  + LioOV 

[Annealed  L2  center]  + Li+  -»  Lin+gOV  . 

They  also  proposed  that  a lithium  ion  can  complex  with  an  annealed 
A-center  to  neutralize  an  L2  center: 

LiOV'[L2-center ] + Li+  ->  ne-  + Lig0V[ Annealed  L2  center  1 

+ Li  — * Lin+20V  . 
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Some  of  these  clusters  are  paramagnetic.  The  variety  of  clusters 
implies  changes  in  electron  binding  energies  and  allows  for  changes  in 
contributions  to  l/Tx  through  modifications  in  the  power  spectra  of  the 
various  donor  electrons  as  well  as  their  concentrations. 

Electron  spin  resonance  and  NMR  results  obtained  for  sample  10-17 
are  in  qualitative  agreement  with  the  cluster  formation  model. 

Figure  17  demonstrates  the  loss  in  the  integrated  electronic  paramag- 
netism with  annealing  and  the  formation  of  diamagnetic  centers.  The 
accompanying  effects  on  the  nuclear  paramagnetism  are  seen  in  Fig.  20. 

A possible  breakup  of  lithium  pairs  (or  clusters)  to  form  isolated 
donor  centers  or  other  paramagnetic  centers  (during  irradiation)  is 
indicated  by  the  suppression  of  the  equilibrium  magnetization  below 
20  K.  During  annealing  the  equilibrium  magnetization  is  observed  to 
recover  towards  its  pre-irradiation  value,  indicating  a possible  re- 
pairing or  re-clustering  of  the  lithium  atoms. 

Measurements  of  Ti  above  room  temperature  for  samples  10-17  and 
0.1-15  after  annealing  indicate  that  the  rate  l/Tx  continues  to  follow 
the  conduction  electron  concentration  and  that  the  variation  of  the 
rate  with  NMR  frequency  is  not  significant.  The  signal-to-noise  ratio 
at  these  temperatures  is  poor  and  Mg  is  not  much  above  background. 

SUMMARY  - EXPERIMENTAL  RESULTS  AND  PHYSICAL  MODELS  - NMR 

In  measuring  the  29Si  spin-lattice  relaxation  time,  T1}  a saturat- 
ing operation  is  first  performed  on  the  sample  in  order  to  reduce  the 
magnetization  M,  due  to  one  29Si  nuclear  moments,  to  zero;  i.e.,  the 
spin  populations  are  equal  for  the  m = \ and  m = - t spin  states.  Then 
the  NMR  field  and  sample  temperature  are  held  constant  for  a period  of 
time  t during  which  the  magnetization  M is  found  to  grow  exponentially 
according  to  the  relation  M(t)  - Mg[l  - exp(-t/T1)].  The  state  of  M(t) 
is  tested  at  time  t by  means  of  a quasi -adiabatic  fast-passage  through 
the  NMR  resonance,  sensing  the  dispersion  mode.  The  process  is  repeated 
for  several  values  of  t.  The  quantity  1 'Ti  then  represents  the  rate  of 
recovery  or  growth  an  M represents  the  steady-state  or  maximum  value 
. 


In  this  investigation  the  rate  of  recovery  l/Tx  is  controlled  en- 
tirely by  either  (l)  conduction  electrons  or  (2)  paramagnetic  electrons 
such  as  single  electrons  occupying  a shallow  or  deep  trap.  There  are 
two  main  factors  involved  in  the  influence  of  the  electrons  on  l/Tx. 

One  of  these,  the  strength  of  the  electron-nuclear  interaction,  can  be 
thought  of  as  the  strength  of  the  local  magnetic  field,  due  to  the 
electron,  at  the  site  of  a Si  nucleus.  This  field  is  called  hj^t). 

Of  course,  different  clei  will  "see"  different  values  so  that  hL  repre- 
sents a mean  of  a stat  stical  variable  with  some  distribution  about  the 
mean.  The  second  factor  in  the  influence  of  the  electron  on  1 /Ti  is 
the  time  variation  of  hL(t)  or  specifically  the  component  in  the  fre- 
quency spectrum  (power  spectrum)  of  hj^t)  which  falls  at  the  NMR  oper- 
ating frequency,  « f. 
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The  power  spectrum  of  the  electronic  local  field  is  related  to  the 
auto-correlation  function  of  hL(t)  through  the  quantity  t called  the 
auto-correlation  time  (or  correlation  time).  A number  of  theories  have 
been  treated  relating  to  this  problem,  of  which  the  most  reasonable 
ones  find  that  the  power  spectrum  and  therefore  l/TL  are  proportional 
to  hLr(l  + of  t ).  The  quantity  t can  then  be  interpreted  as  a measure 
of  the  lifetime  of  an  electron  in  a localized  state  or  in  a spin  state. 
The  frequency  w is  then  taken  to  be  the  NMR  operating  frequency. 

The  above  relation  appears  in  two  different  temperature  ranges  in 
the  analysis  of  the  experimental  results.  In  the  range  from  310  K down 
to  about  60  K,  t is  interpreted  as  the  lifetime  of  an  electron  in  a 
trap.  At  temperatures  below  about  40K,  t is  related  to  the  relaxation 
time  of  the  electron;  i.e.,  the  lifetime  in  a spin  state. 

Upper  Temperature  Region 

In  the  upper  temperature  range  the  lifetime  of  an  electron  in  a 
spin  state,  Te,  is  very  short,  as  is  seen  by  the  very  wide  EPR  linewidth 
observed  at  310  K (Te  less  than  10“‘J  s).  Under  these  conditions  the 
9Si  nucleus  "sees”  the  average  value  of  the  electron  magnetic  moment, 
averaged  over  m = \ and  m = --g  spin  states,  however,  the  average  value 
is  not  zero,  but  rather  q = qe(ueH  /kT)  where  the  fraction  in  parentheses 
comes  from  the  Boltzmann  factor.  The  correlation  time  t in  this  range 
then  represents  the  longer  time  that  an  electron  sits  in  a trap  and 
appears  (to  the  29Si)  to  have  a moment  q (less  than  qe)  and  a local  field 
hL  proportional  to  q.  In  the  plots  of  l/Ti  vs.  T (see  Figs.  5-19), 
samples  0.1-15  and  10-17,  dominated  by  phosphorus  and  lithium,  respec- 
tively, the  measurements  are  seen  to  be  frequency  independent  (converge) 
at  room  temperature,  to  drop  to  lower  rates  as  the  temperature  is 
lowered,  and  to  become  increasingly  frequency  dependent.  This  is  an 
indication  that  un  becomes  increasingly  greater  than  unity  or  t > l/w. 

For  the  UMR  frequencies  of  6,  8.5,  and  10.5  MHz  l/w  has  the  values  of 
28,  20,  and  l6  ns,  respectively,  thus  providing  lower  limits  to  the 
localized  electron  lifetimes.  Below  about  60  K a different  mechanism 
becomes  dominant  as  discussed  below. 

A comparison  of  the  plots  for  0.1-15  and  10-17  shows  that  the  10-17 
plot  has  a smaller  slope  (in  the  60  to  310  K region),  is  shifted  slightly 
toward  lower  temperatures,  and  shows  a slight  frequency  dependence  at 
room  temperature.  This  is  interpreted  to  be  the  result  of  a smaller 
average  trap  binding  energy  for  Li  (lower  temperature)  and  a broader 
distribution  of  trap  depths  about  the  average  (smaller  slope)  including 
some  larger  trap  energies  (room  temperature  frequency  dependence),  thus 
supporting  the  contention  that  Li  in  oxygen-rich  silicon  forms  several 
different  types  of  donor  sites. 

In  the  samples  of  lower  donor  concentrations  (10-15,  1-15,  1-0) 
the  room  temperature  frequency  dependence  is  proportionally  larger  as 
compared  with  the  overall  smaller  rate  l/Ti  (longer  Ti ) indicating  that 
deeper  traps  than  phosphorus  are  competing  on  a roughly  equal  basis  in 
contributing  to  29Si  relaxation.  A projection  of  the  curves  to  higher 
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temperatures  suggests  that  perhaps  600  K would  be  required  to  reach 
frequency  independence,  implying  that  traps  perhaps  twice  as  deep  as 
phosphorus  are  present  in  rougdily  equal  concentrations.  Since  the 
slopes  are  smaller,  a distribution  of  trap  energies  is  indicated. 

The  plot  of  l/Ti  vs.  T for  other  samples,  before  and  after  irra- 
diation sind  annealing,  have  been  studied  in  a similar  manner  and  the 
indications  in  the  upper  temperature  range  (6 0 to  310  K)  are  that,  gen- 
erally, irradiation  causes  a small  decrease  in  temperature-independent 
room  temperature  relaxation  rate,  corresponding  to  a decrease  in  con- 
duction electron  density,  followed  by  a smaller  increase  fter  anneal- 
ing. An  increase  in  frequency  dependence  sind  a smaller  s _ope  is  evident 
alter  irradiation  and  this  effect  does  not  recover  after  annealing. 

This  indicates  that,  generally,  there  is  a broader  distribution  of  trap 
energies  after  irradiation  which  persists  after  annealing. 

Low  Temperature  Region 

In  the  lower  temperature  region  (20  to  40  K)  there  is  a sharp  in- 
crease in  rate  l/Ti , with  decreasing  temperature,  observed  in  all  samples. 

Since  the  electrons  have,  at  this  temperature,  long  lifetimes  as  local- 
ized electrons,  the  power  spectrum  in  this  range  results  from  the  elec- 
tron changing  spin  states  or  from  ground  state  to  excited  state  in  the 
trap  potential.  A number  of  investigators  have  studied  the  electron 
relaxation  rate  1/Te  in  this  region  and  have  found  that  different  pro- 
cesses have  temperature  dependence  proportional  to  T or  T or  exp  A kT. 

These  studies  all  indicate  that  there  is  a very  rapid  increase  in  Te 
with  decreasing  temperature,  and,  if  this  is  taken  to  be  the  correlation 
time  x in  the  expression  h^r/Cl  + w r ) for  the  power  spectrum  of  the 
local  field,  a useful  model  results  for  the  interpretation  of  1 T vs.  T 
in  this  region.  The  position  of  the  characteristic  slope  in  this  region 
is  shifted  toward  lower  temperatures  in  the  case  of  smaller  trap  energies 
(as  in  10-17)  and  shows  a smaller  slope  in  the  case  of  a distribution  of 
trap  energies.  A study  of  the  plots  of  1 'Tj_  vs.  T in  this  range  again 
shows  an  increase  in  average  trap  energy  and  a broader  distribution  of 
energies  after  irradiation  which  does  not  recover  completely  after  anneal- 
ing. 

NM R Magnetization  vs.  Temperature 

Elementary  theory  indicates  that  the  NMR  magnetization  (steady 
state),  Ms,  should  increase  with  decreasing  temperature  in  proportion 
to  1 'T  (Curie's  Law).  The  plot  of  Ms  vs.  T (Figs.  4,  20)  shows  Curie's 
Law  behavior  down  to  temperatures  of  approximately  40  K below  which  a 
sharp  drop  in  Ms  occurs , with  different  samples  dropping  off  at  different 
temperatures.  A comparison  of  the  10-15,  1-15 , and  0.1-15  reveals  that 
higher  P concentrations  da-op  at  higher  temperatuo-es.  In  contrast,  a 
comparison  of  10-17  vs.  (10-15,  1-15,  0.1-15)  vs.  1-0  reveals  that  higher 
Li  concentrations  drop  off  at  lower  temperatures.  This  indicates  (as 
expected)  that  larger  phosphorus  concentrations  correspond  to  higher  N 

concentrations  of  (unpaired)  paramagnetic  electrons  at  low  temperatures 
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which  serve  to  remove  contributions  to  Mo  from  the  observed  WMR  line 
when  some  of  the  electrons  have  t (or  Te)  sufficiently  long  (in  the 
millisecond  range)  to  produce  locally  inhomogeneous  fields. 


The  surprising  result  in  the  contrasting  behavior  of  Ms  vs.  Li 
concentration  is  taken  as  evidence  that  Li,  being  mobile,  tends  to 
form  Li2  pairs  (and  also  Li  P pairs)  or  clusters  which,  at  low  temper- 
atures, become  partially  diamagnetic  thus  reducing  the  low  temperature 
paramagnetic  concentration  with  increasing  Li  concentration.  These 
conclusions  are  supported  by  evidence  obtained  from  EPR  results  and 
from  Hall  effect  and  mobility  studies.  This  effect  is  particularly 
noticeable  in  the  10-17  samples  where  Ms  (with  decreasing  temperature) 
starts  to  decrease  and  then  levels  off  and  starts  to  increase  at  10  K, 
indicating  that  many  of  the  paramagnetic  sites  at  20  K have  become 
diamagnetic  at  10  K. 
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EPR  INVESTIGATIONS 


GENERAL 

Studies  of  the  electron  paramagnetic  resonance  (EPR;  also  called 
electron  spin  resonance,  ESR)  spectra  were  made  to  obtain  supporting 
data  for  the  NMR  study  of  lithium  annealing.  The  main  features  of  the 
EPR  behavior  measured  were  the  spectral  position  or  g-factor,  integrated 
line  amplitude  (magnetization),  linewidth,  and  multiplet  splitting. 

These  quantities  were  obtained  as  functions  of  the  donor  concentrations , 
irradiation  and  annealing  operations,  and  sample  temperature.  The 
samples  were  selected  and  treated  in  essentially  the  same  manner  as  the 
NMR  samples  in  terms  of  irradiation,  annealing,  and  temperature  vari- 
ation. 

Based  on  earlier  reported  studies  the  g-factors  and  multiplet 
splittings  were  used  to  help  identify  the  major  contributing  paramag- 
netic concentrations  affecting  the  NMR  behavior.  The  relative  ampli- 
tudes of  these  spectral  lines,  which  varied  widely  with  temperature, 
donor  concentration,  and  sample  state,  served  to  estimate  the  varying 
concentration  of  NMR  relaxing  centers  (localized  electrons).  The  line- 
width  served  as  an  estimate  of  the  electron  spin-spin  relaxation  rate, 
l/T2e>  which  is  closely  related  to  the  inverse  of  the  electron  spin 
correlation  time,  t,  instrumental  in  the  relaxation  of  the  (Si  nucleus 
at  low  temperatures. 

Two  of  the  most  interesting  developments  resulting  from  the  EHI 
study  are  (l)  the  unusually  large  increase  in  observed  amplitude 
(magnetization)  near  25  K in  the  10-17  lithium-rich  sample  after  irra- 
diation and  (2)  the  consistent  decrease  in  observed  magnetization  below 
20-30  K in  all  samples . 


REVIEW  OF  THEORY 

Only  selected  details  of  the  very  extensive  EPR  literature  and 
theory  will  be  treated  here;  specifically,  those  which  have  a direct 
bearing  on  the  interpretation  of  the  data.  Generally,  the  basic  theo- 
retical developments  are  almost  identical  with  the  basic  theory  for  NMR 
(see  Review  of  Theory  in  preceding  section),  however,  there  are  tradi- 
tionally some  differences  in  notation  and  concepts . 

If  the  electron  moment  interacts  only  with  the  external  field, 
the  spin  Hamiltonian  is 


X g3H-S  (88) 

and  the  resulting  energy  levels  are  given  by 

% g3HM  , (89) 
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where  g is  the  effective  g-factor,  p is  the  Bohr  magneton,  H is  the 
magnitude  of  the  external  magnetic  field,  and  M is  the  magnetic  quantum 
number  and  is  + l/2  or  - l/2  for  the  electron.  A transition  between 
the  levels  requires  a change  in  energy 

AE  - hv  = gpH  (90) 


or 


g - hi  /m  , 


(91) 


where  h is  Planck's  constant  and  v - co/2it  is  the  frequency  corresponding 
to  the  energy  difference.  In  EPR,  v and  H are  measured  at  resonance 
(applied  microwave  frequency  and  external  field)  and  g is  calculated  to 
represent  the  spectral  position,  a characteristic  of  the  specific  spec- 
trum independent  of  the  operating  frequency  and  the  field. 


In  the  event  that  there  are  a number  of  weakly  interacting  spins 
(e.g.,  the  dilute  paramagnetic  system  in  the  present  study)  the  theory 
due  to  Bloch  applies  also  to  the  electron  system  (see  Review  of  Theory 
in  preceding  section).  The  concepts  of  spin-lattice  and  spin-spin  re- 
laxation times  for  the  electron,  TLe  and  T2e,  are  carried  over  without 
change  from  the  NMR  theory  although  their  range  of  values  is  normally 
smaller . The  theoretical  machinery  for  dealing  with  the  electronic 
magnetization  in  the  rotating  coordinate  system  in  terms  of  the  effec- 
tive magnetic  field  is  the  same. 

If  the  electron  also  interacts  (hyperfine  interaction,  scalar  form) 
with  a nuclear  moment,  such  as  the  localized  electron  and  the  phosphorus 
nucleus,  the  spin  Hamiltonian  is 


gPH-S  + AI-S  , (92) 

where  I is  the  nuclear  spin  and  A is  the  strength  of  the  electron- 
nuclear  scalar  interaction,  proportional  to  the  electron  wave -funct ion 
density  at  the  nucleus.  The  energy  levels  (first  order)  are  given  by 


%,m  “ gP™  + AMn  • 

A transition  involving  ^ + 1 requires  energy  changes  of 


(93) 


which,  in  the  case 


AEj,! 

of  phosphorus 
hv+ 


gPH  + mA  , 

(spin  = l/2),  becomes 
- gPH  ± A/2 


(94) 

(95) 


and  two  resonances  are  observed, 
in  gauss,  provided  A/h  is  larger 
be  resolved. 


separated  by  A/h  in  frequency  or  A/’gP 
than  the  linewidth  so  that  they  can 
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If  two  phosphorus  sites  are  sufficiently  near  one  another  (less 
than  ~ 50  A or  10  unit  cell  dimensions  in  silicon)  at  low  temperature, 
the  two  localized  electrons  interact  with  one  another  (exchange  inter- 
action's well  as  with  the  phosphorus  nuclei.  The  spin  Hamiltonian  is 

"K  = gp(SL  + S2)-H  + JSL  -S2  + A(Ll  + IaMSj.  + S2)  (96) 

for  A « J,  the  electrons  form  a resultant  spin  §,  - Si  + S2  with  S = 0 
or  1,  and  the  electron  pair  forms  a singlet  state  (S  0)  or  a triplet 
state  (S  l) . The  singlet  state  is  nonpar amagnetic  and  does  not  con- 
tribute to  the  EPR  signal  nor  to  NMR  or  EPR  relaxation;  i.e.,  the  forma- 
tion of  singlet  states  tends  partially  to  "clean  up"  the  paramagnetic 
impurities  in  a sample.  This  was  observed  to  occur  in  this  study  (see 
Fig.  21)  in  the  25  K range  and  below.  The  resulting  energy  levels  for 
the  triplet  state  are  (first  order) 

%,m  = 80HM  + i AMm  + J J , (97) 

where 

m = mjL  + m2  . 


For  AM  ± 1,  transition  energies  are 

= g@H  + | Am  . 

Possible  values  for  m - mx  + m2,  for  two  phosphorus  nuclei,  are  1,  0, 

-1  with  the  probability  of  m = 0 twice  as  large  as  m - 1 or  m = -1. 

The  resulting  spectrum  consists  of  two  outer  lines  separated  by  A with 
a third  line  with  double  amplitude  centered  between  the  outer  lines . 

The  center  line  has  the  same  g-factor  as  the  unsplit  line. 

If  there  are  three  interacting  electron- phosphorus  sites,  the 
spectrum  consists  of  four  equally  spaced  lines  centered  as  before,  with 
the  outer  two  coinciding  with  the  one-electron  doublet  and  with  the 
two-electron  outer  pair.  The  amplitudes  have  ratios  of  1:3: 3=1  in  the 
three-electron  system.  The  pattern  extends  to  higher  numbers  of  elec- 
trons . Figure  21  shows  a superposition  of  a strong  one-electron 
doublet,  a weaker  two-electron  triplet,  and  an  even  weaker  three- 
electron  quadruplet. 

The  energy  level  for  the  singlet  state,  S = 0,  M 0,  is 

E - 3J/4  , (98) 

which  is  lower,  by  the  amount  J,  than  the  average  of  the  triplet  state 
levels.  For  phosphorus  donor  sites  separated  by  ~ 25  A or  less  the 
value  of  J is  in  the  same  range  of  values  as  the  electron  binding 
energy  to  an  isolated  phosphorus  donor  site  so  that  the  actual  trap 
levels  in  donor  pairs  are  shifted  significantly.  At  room  temperature 
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the  levels  are  almost  equally  likely  to  be  occupied  so  that  the  average 
trap  energy  is  almost  unchanged  from  a single  donor  energy.  At  lower 
temperatures,  as  kT  becomes  less  than  J,  the  lower  level  nonparamsgnetic 
singlet  tends  to  have  a higher  population  resulting  in  a decrease  in 
paramagnetism. 

In  the  case  of  phosphorus,  the  internuclear  distances  are  fixed 
(phosphorus  is  not  very  mobile),  however,  if  the  two  (or  more)  donors 
are  mobile  at  room  temperature,  such  as  lithium,  the  corresponding 
exchange  interaction  between  two  Li,  Li  and  LiO,  or  between  Li  and  P 
serves  to  provide  one  more  mechanism  for  tending  to  bind  mobile 
impurities  to  specific  sites.  Since  the  binding  is  weak  it  should  be 
easy  to  break. 

The  Li  nuclear  hyperfine  interaction  with  the  localized  electron 
(92.6%  abundant  Li  spin  3/2;  7.4%  abundant  Li  spin  - l)  is  expected 
in  theory  to  be  similar  to  that  of  phosphorus,  however,  the  interstitial, 
lower  symmetry,  position  of  Li  leads  to  a greater  p-orbital  character 
and  a smaller  electron  density  at  the  nucleus,  resulting  in  smaller 
splitting  factor,  A.  Furthermore,  Li  is  expected  to  be  found  in  a 
variety  of  bonding  situations  so  that  the  splitting  is  unresolved  and 
the  Li-related  EPR  line  is  somewhat  broadened. 


EXPERIMENTAL  PROCEDURE 


■ 

ft* 

*•4 

[ 

*8 

h 

*. . 

4 

* 

b 

rW 

K 

s 

7 

u 

i 

i 

J 

The  experimental  arrangement  for  the  EPR  study  is  very  similar  to 
commercial  systems  and  is  essentially  the  same  as  the  NMR  system  (see 
Experimental  Procedure  in  preceding  section),  except  that  a klystron, 
wave  guide,  and  microwave  cavity  replaces  the  transmitter,  and  the  rf 
signal  is  amplitude  detected  without  rf  amplification.  The  field  is 
modulated  at  100  kHz  and  the  resulting  signal  at  this  frequency  is 
amplified  and  phase-detected  to  provide  the  derivative  of  the  EPR 
absorption  line  for  recording.  The  low-temperature  and  temperature- 
control  systems  are  essentially  the  same  as  for  NMR. 

The  samples  for  the  EPR  study  were  cut  from  one  of  the  NMR  slabs 
to  4 x 10  x 0.25  mm,  a volume  of  0.01  cm',  so  that  the  number  of  elec- 
tron spins  in  the  sample  was  10!L4-1015.  The  sample  was  held  in  the 
Dewar  by  means  of  a 4 mm  quartz  tube,  slotted  to  contain  the  slab  with 
the  long  dimension  vertical.  The  sample  could  be  rotated  about  the 
vertical  with  respect  to  the  horizontal  external  field.  An  auxiliary 
sample  of  DPPH  was  located  in  the  cavity,  outside  the  Dewar,  so  that  it 
could  easily  be  removed  or  re-inserted,  to  provide  a reference  marker 
both  for  g-factor  and  for  amplitude;  thus  calibrated  amplitude  data 
were  available.  Since  the  DPPH  was  always  at  room  temperature  the  cali- 
brating reference  was  not  affected  by  changes  in  sample  temperature. 

The  investigative  procedure  consisted  of  studying  the  selected 
samples  over  a range  of  temperature  from  6 to  K,  recording  the 
observed  spectra  along  with  the  reference  spectrum,  sample  temperature, 
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field  markers,  and  other  pertinent  data.  In  certain  cases,  the  sample 
was  rotated  to  display  possible  orientation-dependent  spectra  but  none 
were  noted.  The  derivative  peak-to-peak  amplitude  and  linewidth  were 
measured  ar i used  to  recover  absorption  line  area  data.  The  absorption 
line  area  is  proportional  to  the  product  of  the  amplitude  and  the  square 
of  the  width  of  the  derivative  line  with  the  proportionality  constant, 
dependent  upon  the  line  shape,  Gaussian  or  Lorentzian.  The  variation 
of  this  integrated  amplitude,  proportional  to  rf  magnetization,  Mg,  was 
then  studied  as  a function  of  sample  temperature.  The  g- factor  and 
splitting  factors  were  noted  as  were  the  critical  temperatures  at  which 
changes  in  line  structure  occurred,  such  as  single  line  to  doublet. 


EXPERIMENTAL  RESULTS  AND  ANALYSES 

The  sample  magnetization,  Me,  in  all  cases  showed  an  increase  after 
irradiation  followed  by  a partial  recovery  after  annealing.  The  greatest 
increase  occurred  in  the  10-17  sample  (Fig.  17)  in  the  25  K temperature 
region. 

The  absorption  lines  do  not  become  very  strong  and  narrow  until 
the  temperature  is  lowered  to  ~ 100  K.  This  indicates  that  the  spin- 
spin  relaxation  time,  Taes  is  less  than  about  2 ns  in  the  upper 
temperature  region,  and  thus  supports  the  conclusion  that  the  frequency- 
dependent  NMR  relaxation  in  that  region  is  due  to  the  longer  localized- 
electron  residence  time  (trap  lifetime)  and  also  that  the  electronic 
moment  (q)  average  is  equal  to  qe(peH/kT),  as  in  Eq.  (4-9),  for  times 
comparable  with  the  nuclear  spin-spin  relaxation  time  Ta. 

As  the  temperature  is  decreased,  the  magnetization,  Mg,  shows 
generally  an  increase  related  to  the  Cur^e-law  modified  by  the  fact 
that  the  number,  Np,  of  paramagnetic  electrons  is  also  increasing,  in 
addition  to  the  Boltzmann  factor  (qeIl/kT)  responsible  for  the  Curie-law. 
In  the  case  of  the  phosphorus  donor  electron  EPR  line,  a variation  from 
the  Curie-law  behavior  occurs  in  the  range  near  35  K because  of  the 
change  from  a single  line  above  40  K to  a doublet  at  30  K.  This  be- 
havior allows  an  estimate  of  lifetimes  as  follows:  at  about  40  K where 

the  doublet  lines  begin  to  develop,  21  G or  370  x 10'  rad/s  away  frcm 
the  singlet,  the  electron  lifetime  in  a localized  state  at  the  phosphorus 
site  is  ~ 1/370  qs.  Then  at  ~ 30  K where  the  individual  line  in  the 
doublet  narrows  to  ~ 30,  the  lifetime  is  longer  by  a factor  of  7 or 
~ 1/53  qs.  These  values  for  lifetime  of  states  are  mean  values  of 
statistical  distributions  which  are  as  broad  as  the  relative  variation 
in  electron-electron  interaction  energies  in  the  lattice  so  that  sub- 
stantial numbers  of  electrons  have  longer  (or  shorter)  corresponding 
lifetimes  by  factors  of  ~ 10  (or  1/10). 

At  temperatures  in  the  20-25  K range  the  electron  magnetization, 

Mg,  consistently  begins  to  decrease  with  decreasing  temperature,  con- 
trary to  Curie-law  behavior.  This  appears  to  result  from  at  least  two 
causes;  i.e.,  (l)  pairs  and  clusters  of  donor-site  electrons  interact 
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and  become  partially  diamagnetic  (nonparamagnetic)  and  (2)  long  lifetime 
electron  states  (longer  than  the  reciprocal  linewidth)  increase  internal 
field  inhomogeneities  and  remove  contributions  from  the  rf  magnetization 
as  observed  by  EPR. 

The  reduction  in  magnetization  below  ~ 25  K occurs  generally  at 
lower  temperatures  in  lithium-rich  samples  indicating  that  lithium- 
related  trap  energies  are  smaller  than  phosphorus . The  ratio  appears 
to  be  about 

In  the  lithium-rich  10-17  sample,  the  electronic  magnetization, 

Me,  was  observed  to  increase  by  an  unusually  large  amount  after  irra- 
diation (before  annealing),  especially  in  the  temperature  range  near 
25  K (Fig.  17). 

It  is  believed  that  this  Me  increase,  by  irradiation,  over  the 
unirradiated  state  may  result  partially  from  breaking  up  lithium- 
related  pairs  and  clusters  which  would  normally  be  partially  diamagnetic, 
thus  increasing  the  paramagnetism,  especially  at  lower  temperatures 
where  shallow-trap  electrons  become  strong  Mg  contributors.  Of  course, 
it  is  expected  that  deeper-level  paramagnetic  electrons  are  also  created 
by  bombardment  but  these  electrons  are  Mg  contributors  at  higher  tem- 
peratures, even  at  room  temperature. 


SUMMARY 

At  room  temperature  the  EPR  line  is  very  weak  and  broad,  if  seen 
at  all,  indicating  a very  short  lifetime  of  the  electron  in  a paramag- 
netic state  (one  ns  or  less).  At  about  100  K the  line  has  narrowed  to 
about  3 G (Tog  about  10  ns)  and  the  g-factor  can  be  measured  to  better 
than  one  part  in  10’ . In  P-rich  samples,  g corresponds  to  reported 
values  for  phosphorus  site  electrons.  In  Li-rich  samples,  two  different 
g factors  were  observed  and  were  in  agreement  with  g for  LiO  and  LiOV 
type  sites. 

At  temperatures  lower  than  30  K the  P-site  line  was  observed  to 
split  into  a doublet  (aG  = 42  G)  and  in  P-rich  samples  weaker  triplet 
and  quadruplet  structure  was  observed.  This  splitting  behavior  results 
from  interaction  of  the  electron  with  the  phosphorus  nucleus  (s  = |)  in 
the  case  of  the  42  G doublet.  The  triplet  and  quadruplet  result  from 
an  exchange  interaction  between  two  or  three  adjacent  phosphorus  sites 
and  appear  when  kT  has  become  smaller  than  the  exchange  energy  J.  The 
two  electrons  can  assume  a triplet  state  (S  = l)  and  remain  paramagnetic, 
or  a singlet  state  (S  = 0)  and  become  diamagnetic.  In  the  latter  case 
the  contribution  is  lost  from  EPR  magnetization  and  also  NMR  magnetiza- 
tion is  increased  and  NMR  relaxation  rate  is  decreased.  In  the  triplet 
state  the  two  phosphorus  nuclear  moments  combine  to  produce  the  three- 
component  line  in  the  EPR  spectrum.  The  lithium  lines  do  not  show  a 
nuclear  hyperfine  splitting  even  though  the  dominant  Li  nucleus  has  a 
spin  of  3/2.  This  results  from  the  lower  symmetry  of  the  interstitial 
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Li  position  and  the  strong  p-state  character  of  the  donor  electron  as 
contrasted  with  the  strong  s-state  character  of  the  phosphorus  donor 
electron.  In  the  p-state  the  electron  density  at  the  nucleus  is  small 
and  thus  also  the  hyper fine  interaction. 

The  pairing  of  Li  (or  clustering)  is  revealed,  however,  in  the 
behavior  of  the  10-17  sample  upon  irradiation.  An  increase  in  electron 
magnetization  of  an  order  of  magnitude  was  observed  at  25  K after  irra- 
diation, which  then  reduced  nearly  to  the  original  level  after  24  hr 
room  temperature  annealing.  This  can  be  accounted  for  by  a break-up  of 
lithium  clusters  by  the  irradiation,  followed  by  a relatively  rapid 
return  to  cluster  formation  upon  annealing.  This  model  is  also  sup- 
ported by  the  Hall  measurements  and  by  the  NMR  Ms  vs.  T data. 

In  contrast  with  Curie  Law  behavior,  all  of  the  samples  generally 
show  a decrease  in  electronic  steady-state  rf  magnetization,  Me,  with 
decreasing  temperatures  in  the  region  below  about  25  K.  Several  mecha- 
nisms can  contribute  to  this  effect:  (l.)  change  in  neighboring  pairs 

from  paramagnetic  to  diamagnetic  states  when  kT  becomes  less  than  the 
pair  exchange  energy  J,  (2)  increase  in  local  field  inhomogeneity  when 
the  electronic  relaxation  time  becomes  sufficiently  long,  or  (3)  partial 
saturation  of  the  spin  system.  Tests  were  made  to  see  if  partial  satu- 
ration could  be  present  but  the  indications  were  that  this  was  not 
responsible  for  the  observed  non-Curie  behavior.  Evidence  from  EPR, 

NMR,  and  Hall  measurements  all  tend  to  support  the  first  two  of  the 
above  mechanisms.  Thus,  if  two  (or  more)  donor  sites  are  neighbors 
(30-40  A or  less),  whether  they  are  P-P,  Li -Li,  P-Li,  Li-LiOV,  or  other, 
at  low  temperatures  the  electrons  localized  on  these  sites  form  exchange- 
energy  pairs,  some  of  which  are  in  the  diamagnetic  lowest  energy  state 
and  do  not  contribute  to  the  EPR.  magnetization  Me.  In  addition,  the 
remaining  paramagnetic  electrons  achieve  long  relaxation  times  suffi- 
cient to  increase  the  local  field  inhomogeneity  and  remove  contributions 
from  the  observed  EPR  linewidth. 

The  increase  in  low-temperature  EPR  magnetization.  Mg,  after 
electron  irradiation  is  interpreted  to  indicate  that  these  pairs  (or 
clusters)  are  scattered  or  broken-up  by  the  bombardment,  and  further 
that  the  clusters  tend  to  reform  after  annealing,  again  reducing  Mg. 


HALL  EFFECT  AND  RESISTIVITY 


GENERAL 


For  verification  of  resistivity  and  conduction  electron  concentra 
tion,  Hall  coefficients,  and  resistance,  measurements  were  made  on  the 
selected  samples  under  varying  conditions  of  sample  temperature,  irra- 
diation, and  annealing. 


Before  irradiation  the  measurements  were  in  reasonably  good  agree- 
ment with  the  stated  resistivities  and  Li  doping  levels.  The  behavior 
after  irradiation  was  about  as  expected  in  all  except  the  10-17  lithium 
rich  sample  which  showed  an  increase  in  conduction  electron  concentra- 
tion upon  irradiation. 


EXPERIMENTAL  PROCEDURE 


Both  for  Hall  voltage  and  resistance  measurements,  a constant- 
current  dc  source  was  devised  that  held  the  Hall  current  constant  to 
better  than  0.5 % over  sample  resistance  changes  (with  varying  tempera- 
ture) of  more  than  three  orders  of  magnitude.  Available  constant 
currents  ranged  from  0.1  uA  to  10  mA.  The  Hall  voltage  and  the  resis- 
tance were  then  monitored  by  accurate  digital  voltmeters  with  high  input 
impedances . Both  Hall  current  and  external  field  directions  were 
reversed  independently  in  order  to  minimize  spurious  effects.  Consist- 
ent readings  were  obtained  for  conduction  electron  concentrations  down 
to  1012  cm-3  and  below. 


Sample  slabs  20  x 10  x 0.25  mm  were  etched  and  provided  with  five 
gallium-indium  eutectic  alloy  contact  pads,  with  the  use  of  ultrasonic 
vibration  at  room  temperature.  One  contact  was  located  at  each  end,  a 
third  centered  on  one  edge,  and  two  others  were  on  the  opposite  edge 
one-third  the  distance  from  each  end.  Contact  with  these  pads  was  made 
by  spring-loaded  pressure  to  brass  contacts  on  the  holder  and  through  a 
cable  to  the  meters.  The  edge  contacts  were  used  to  develop  a Hall 
voltage,  balanced  against  zero-field  voltage  drops,  or  the  two  on  a 
common  edge  were  used  to  develop  a resistance-related  voltage.  The 
end  contacts  were  used  to  apply  the  Hall  current  and  also  to  measure 
two-terminal  resistance.  The  latter  technique  was  found  to  be  satis- 
factory for  most  measurements  in  this  study. 


The  magnetic  field  was  provided  by  the  NMR  electromagnet  at  10  kG. 
Also,  the  same  NMR  Dewar  and  control  system  was  used  to  vary  the  sample 
temperature.  After  irradiation  at  room  temperature,  the  samples  were 
stored  at  77  K until  measurements  were  started.  The  sample  was  then 
transferred  to  the  test  Dewar  (~  15  min)  and  low-temperature  measure- 
ments were  made  first.  The  samples  were  again  examined  after  annealing 
~ 60  hours . 
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EXPERIMENTAL  RESULTS  AND  ANALYSES 

The  resistivity  and  conduction  electron  concentration,  observed  in 
the  various  samples  before  irradiation,  agree  very  well  with  values 
expected  on  the  basis  of  the  stated  phosphorus  and  lithium  concentra- 
tions. These  quantities  also  follow  reasonably  well  the  changes 
expected  with  temperature  variation.  In  addition,  the  expected  change 
to  lower  conduction  electron  concentration  after  irradiation  was  found 
in  all  except  the  10-17  lithium-rich  sample  above  100  K.  The  expected 
change  to  higher  conduction  electron  concentration,  n,  after  annealing 
was  found  in  all  lithium-containing  samples,  including  the  10-17 • 

Figure  22  shows  the  changes  found  in  n in  the  0.1-15  phosphorus- 
rich  sample.  A decrease  after  irradiation  at  room  temperature  ~ 25% 
increases  to  larger  percentage  at  low  temperatures,  indicating  an  in- 
crease in  average  trap  depth  below  the  conduction  band  (electrons 
become  localized  at  higher  temperatures). 

The  plot  of  mobility  vs.  temperature  for  sample  0.1-15  (Fig.  23) 
shows  an  unusual  behavior.  After  irradiation  the  mobility  is  lower  in 
the  middle  temperature  range  and  higher  in  the  low  range,  than  either 
before  irradiation  or  after  annealing,  the  cross-over  being  in  the  70  K 
region . This  behavior  was  found  to  occur  in  all  samples . It  appears 
that  a wider  variety  of  trap  energies  is  present  after  irradiation  than 
before  irradiation  or  after  annealing,  including  some  very  shallow 
levels . 

The  10-17  conduction  electron  concentration  (Fig.  24)  in  the  upper 
temperature  range,  shows  the  unusual  behavior  of  increasing  after  irra- 
diation. This  suggests  that  a variety  of  trap  energies  was  produced 
with  more  shallow  than  deep  compared  with  unirradiated  conditions . 
Annealing  the  sample  then  restored  some  of  the  deep  traps  to  shallow 
traps  to  further  increase  the  conduction  electron  concentration. 

The  mobility  of  the  10-17  sample  (Fig.  25)  shows  again  the  unusual 
behavior,  after  irradiation,  found  in  the  other  samples.  It  is  more 
pronounced,  however,  in  the  10-17  sample. 

The  other  samples  of  lower  donor  impurity  concentrations  all  show 
changes  similar  to  the  0.1-15  but  involving  smaller  changes,  less  than 
10fo. 


SUMMARY 

Conduction  electron  density,  n;  resistivity,  p;  and  mobility,  p, 
were  obtained  from  Hall  effect  and  resistance  measurements  as  a function 
of  temperature  and  then  redetermined  after  electron  irradiation  and 
after  annealing  (see  Figures  22-25). 
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Before  Irradiation 
After  Irradiation 
After  Annealing 


TEMPERATURE  (K) 


Figure  22.  Conduction  electron  concentration  for  sample  0.1-15 


Before  Irradiation 
After  Irradiation 
After  Annealing 


TEMPERATURE (K) 


Figure  25.  Mobility  vs.  temperature  for  sample  10-17 


v Generally,  the  conduction  electron  density  showed  a decrease  after 

irradiation,  indicating  a loss  of  conduction  electrons  to  deep  traps, 
except  in  the  case  of  the  10-17  sample  (lithium- rich)  where  the  conduc- 
tion electron  density  showed  an  increase  after  irradiation  as  well  as 
after  annealing.  This  is  viewed  as  supporting  evidence  for  lithium 
pairing  or  clustering.  After  irradiation,  the  deeper  trap  defects  are 
outnumbered  in  the  lithium-rich  sample  by  the  scattered  Li  interstitial 
shallow  donors  thus  increasing  the  conduction  electron  density.  Upon 
annealing,  the  further  increase  in  conduction  electron  density  indicates 
that  the  lithium  heals  the  deep  trap  defects  more  rapidly  than  it  re- 
turns to  clusters. 


The  mobility,  generally,  shows  an  increase  with  decreasing  tempera- 
ture, however  the  greatest  rate  of  increase  occurs  after  irradiation 
(before  annealing),  especially  in  the  lithium-bearing  samples.  Since 
the  low-temperature  mobility  reflects  the  contribution  from  shallow- 
trap  electrons  (the  only  ones  still  ionized),  the  greater  low- temperature 
mobility  after  irradiation  apparently  results  from  very  shallow  traps, 
again  supporting  the  increase  in  isolated  Li  interstitial  donors  upon 
irradiation. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A number  of  the  more  important  conclusions  reached  in  the 
investigation  are  summarized  in  this  section. 

(1)  The  technique  of  utilizing  NMR  as  an  investigative  tool  for 
following  radiation  damage  and  lithium  annealing  effects  in  bulk  n-type 
silicon  has  proved  useful  in  several  ways. 

(2)  The  relaxation  rate  of  the  23Si  nucleus  follows  the  conduction 
electron  concentration  in  the  neighborhood  of  room  temperature. 

(3)  At  somewhat  lower  temperatures , residence  times  of  conduction 
electrons  in  trap  sites  are  revealed  by  the  NMR  frequency  dependence  of 
the  relaxation  rate. 

(4)  Variations  of  electron  trap  energy  levels  from  the  normal 
phosphorus  level  are  deduced  in  high  resistivity  samples  and  in  samples 
after  irradiation  from  reductions  in  the  l/TL  vs.  temperature  slope. 

(5)  The  correlation  time  of  electron  spin  states  and  for  transi- 
tions between  ground  state  and  excited  state  localized  electrons  follow 
from  the  behavior  of  l/T]_  in  the  lower  temperature  range. 

(6)  Variations  from  Curie-law  behavior  are  seen  in  the  measurement 
of  magnetization  vs . temperature  supporting  the  presence  of  donor  site 
pairs  and  clusters . 

(7)  EPR  and  Hall  results  also  support  his  model. 

(8)  The  behavior  of  lithium  under  irradiation  and  annealing  appears 
to  be  related  to  the  presence  of  oxygen  in  these  crucible-grown  samples. 

(9)  The  frequency  dependence  at  room  temperature  in  l/Tx  indicates 
the  presence  of  deeper  than  normal  donor  sites  in  high  resistivity 
samples . 

(10)  Donor  sites  shallower  than  normal  appear  to  be  present  in  the 
10-17  lithium-rich  sample  after  irradiation. 

(11)  The  reduction  of  electron  paramagnetism  below  Curie-law 
predictions  is  shown  and  related  in  part  to  pairing  or  clustering. 

(12)  The  changes  in  1/Ti  vs.  T after  irradiation  and  after  anneal- 
ing indicate  that  deeper  level  donor  sites  are  introduced  and  then 
reduced  in  number  but  not  completely  removed  by  annealing.  In  the 
lithium-rich  sample  it  appears  that  annealing  of  deeper  level  damage 
sites  proceeds  more  rapidly  than  the  return  of  Li  to  pairs  and  clusters. 
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The  most  obvious  recommendation  involves  the  study  of  similarly 
doped  samples  of  oxygen-lean  n-type  silicon  as  well  as  samples  with 
moderate  oxygen  content  in  order  to  gain  more  insight  into  the  role  of 
oxygen  in  lithium  annealing.  A proposal  to  this  effect  has  been  sub- 
mitted. 

Additional  theoretical  work  is  needed  regarding  the  energy  con- 
siderations surrounding  donor  pairs  and  clusters . 
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An  improved  theoretical  model  may  be  possible  for  describing  the 
upper  temperature  frequency  dependence  of  l/TL  resulting  from  short 
lifetime  visits  of  conduction  electrons  to  localized  trap  levels. 

The  behavior  of  “3Si  relaxation  processes  in  p-type  silicon  should 
be  studied  for  purposes  of  comparison  as  well  as  for  adding  to  the 
understanding  of  that  material . 
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Table  3.  Sample  1-15  (Before  Irradiation) 
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Table  5-  Sample  10-17  (Before  Irradiation) 


9m 


Table  9*  Sample  10-17  (After  Irradiation,  Before  Annealing) 
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